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I.  INTRODUCTION 


The  Ballistic  Research  Laboratory  (BRL)  of  the  U.S,  Army  Armament 
Research  and  Development  Command  (ARRADCOM)  has  a  wide  interest  in  the 
shaped-charge  problems,  ranging  from  detailed  studies  of  the  flow  char¬ 
acteristics  of  the  collapsing  liner,  to  designing  practical  devices  for 
future  warhead  applications.  For  these  efforts  several  experimental 
and  theoretical  techniques  are  employed.  Often,  it  is  necessary  to 
determine  parametric  relationships  in  the  design  application.  Since 
it  is  more  feasible  to  utilize  theoretical  calculations  that  are  eco¬ 
nomically  sound  for  parametric  studies,  the  BRL  has  several  finite- 
difference,  hydrodynamic,  computer  codes  that  have  been  applied  to 
shaped-charge  problems.1*2  Although  these  codes  are  adaptable  to  vari¬ 
ous  geometrical  considerations,  they  require  operator  experience  and  a 
seasoned  analyst  to  insure  proper  application.  Further,  large, 
high-speed  computers,  and  long  calculational  times  are  necessary. 

Quite  often,  it  is  desirable  to  have  a  simplified  procedure  for 
addressing  parametric  design  studies  quickly  and  economically.  The 
BRL  computei'  code  named  BASC  (BRL  Analytical  Shaped  Charge) 3  was  formu¬ 
lated  from  analytic  expressions  to  provide  such  capability.  Although 
several  advantages  occur  with  the  original  BASC  approach,  several  areas 
of  difficulty  were  experienced,  particularly  those  relating  to  accurate 
calculation  of  jet  tip  or  lead  pellet  behavior  and  confined  charges. 
Extensive  semi-empirical  functions,  regarding  liner  acceleration  and 
confinement  effects,  have  been  included  to  provide  more  accurate  repre¬ 
sentations.  This  improved,  simplified  procedure,  hereinafter  referred 
to  as  the  BASC  code  also,  together .with  additional  refinements  are  pre¬ 
sented  and  discussed  in  this  report. 


lJ.  T.  Harrison  and  H.  R.  Karpp,  " Terminal  Ballistic  Applications  of 
Hydrodynamic  Computer  Code  Calculations BRLR  19L4,  April  1977. 

JAD  MO 4 106,5)  „  . 

V.  T.  Harrison ,  "A  Comparison  Between  the  Euler  van.  Hydrodynamic 
Computer  Code  ( BRLSC )  and  Experimental  Collapse  of  a  Shaped  Charge 
Liner,"  ARBRL-MR-02841 ,  June  1978  (AD  AO 597 11) 


3J.  Harrison,  R.  DiPevsio ,  R.  Karpp  and  R.  Jameson ,  "A  Simplified 
Shaped  Charge  Computer  Code:  BASC,"  DEA-AF-F/G-7304  Technical 
Meeting:  Physics  of  Explosives,  Vol  II,  April-May  1974,  Paper  13 
presented  at  the  Naval  Ordnance  Laboratory,  Silver  Spring,  MD. 
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The  BASC  code  is  an  assembly  of  various  theoretical  and  empirical 
techniques.  Central  to  the  procedure  utilized  in  BASC  is  the 
Defoumeaux  model4  for  final  plate  velocity  resulting  from  the  shock 
of  an  adjacent  detonating  explosive.  This  assumption  is  adequate  for- 
portions  of  the  liner  which  are  initially  removed  (remote)  from  the 
collapse  (jet  formation)  region  or  cone  axis.  In  actuality,  several 
shock  reverberations  are  required  to  achieve  a  final  liner  metal  velo¬ 
city.  Material  near  to  the  apex  of  the  cone  can  enter  the  collapse 
process  long  before  the  liner  is  accelerated  fully  and,  hence,  does 
not  achieve  its  ultimate  velocity.  This  leads  to  the  well  known  pheno¬ 
mena  referred  to  as  "the  inverse  velocity  gradient"  which  forms  the 
massive  ^et  tip.  This  phenomena  has  been  observed5  and  calculated 
earlier.’ »2«6  The  author  has  modified  the  Defourneaux  model  to  account 
for  the  time-dependent  acceleration  of  the  liner  resulting  in  a  gradual 
build-up  to  the  ultimate  collapse  velocity.  Since  liner  elements  near 
the  apex  region  of  a  cone  will  not  achieve  this  ultimate  velocity,  the 
first  element  of  the  jet  will  move  more  slowly  than  the  following  clc- 
onts  of  the  jet.  The  faster  elements  collide  and  are  compressed  into 
the  massive  lead  pellet.  The  final  jet-tip  velocity  will  become  the 
mass-weighted  average  of  these  inverse  velocity  elements.  This  is  the 
jet-tip  velocity,  which  will  be  used  in  jet  penetration  theory. t  BASC 
uses  a  combination  of  the  shaped-charge  penetration  of  DiPersio,  Simon 
and  Merendino  (DSM)7  and  the  piece  wise  penetration  of  Defourneaux4. 

The  DSM  model  is  used  to  calculate  total  penetration-standoff  curves 
and  the  piecewise  penetration  model  is  used  to  calculate  whole  profiles. 

In  addition  to  the  inverse  velocity  gradient,  flow  into  the  stagna¬ 
tion  region  during  the  collapse  process  from  particles  near  the  apex  of 
a  conical  liner  may  be  supersonic  and  fail  to  form  a  jet  or  form  a  so- 
called  incoherent  jet.  This  is  the  jet-no-jet  criteria.  The  criteria 
used  in  BASC  is  a  modified  version  of  the  supersonic  limitations  of 
Chou,  ct.  al.8  The  final  theoretical  technique  used  in  BASC  is  the 


4M.  Defoumeaux ,  "Hydrodynamic  Theory  of  Shaped  Charges  and  of  Jet 
Penetration Memorial  DeL'art  Ille’rie  Franaasise-T ,  44,  1970. 

SR.  DiPersio ,  C.  W.  Whiteford,  and  J.  Simon,  "An  Experimental  Method 
of  Obtaining  Collapse  Velocities  of  the  Liner  Walls  of  a  Linear 
Shaped-Charge  Liner,"  BRL-MR-1696,  September  1965. (AD  # 478326 ) 

6/l.  Kiwan  and  H.  Wisniewski,  "Theory  and  Computations  of  Collapse  and 
Jet  Velocities  of  Metallic  Shaped  Charge  Liners,"  BRLR-1620 , 

November  1972.  (AD  H907161) 

7R.  DiPersio,  J.  Simon ,  and  A.  Merendino,  "Penetration  of  Shaped- 
Charge  Jets  Into  Metallic  Targets,"  BRLR  1296,  September  1965.  (AD  # 476717 ) 

8P.  Chou,  J.  Carleone,  R.  Karpp,  "Criteria  for  Jet  Formation  from 
Impinging  Shells  and  Plates,"  J.  Appl.  Physics ,  Vol.  47,  No.  7, 

July  1976. 


Pugh,  Eichelberger,  and  Rostoker  (PER)  theory  of  jet  formation.9 

This  report  includes  a  description  of  analytical  equations,  the 
procedures  for  their  application  and  a  FORTRAN  listing  of  BASC.  Ex¬ 
perimental  results  for  selected  shaped  charge  problems  are  presented  in 
a  comparison  to  the  calculations.  The  range  of  useful  applications  is 
discussed  as  well  as  the  limitations  of  the  approach. 

BASC  is  a  "living  code"  which  has  provided  insights  to  the  jet- 
formation  process  and  is  a  very  good  tool  for  parametric  design  for  a 
selected  class  of  problems. 

A  simplified  flow  chart  and  the  FORTRAN  code  listing  is  presented 
in  Appendices  A  and  B,  respectively. 


II.  GOVERNING  EQUATIONS 

The  initial  equation  used  in  the  liner  acceleration  portion  of 
the  BASC  code  determines  the  angle  of  liner  bending,  <t>,  produced  by 
a  detonation  wave  traveling  with  a  velocity,  D,  and  inclined  to  the 
liner  wall  at  an  angle,  i  (the  angle  of  incidence).  This  relationship 
is  illustrated  in  Figure  1.  The  author  has  modified  the  original 
platepush  relationship  of  Defourneaux  to  be 


1  1 
♦  “ 


K  £* 
eB 


B  *  1+A/p  e 
c  ^ 


(1) 


where  p  and  e  are  the  density  and  thickness  of  the  liner  wall,  and 
e  is  the  explosive  thickness.  Added  are  p  ,  oc>  and  A  which  are, 

respectively,  density,  thickness,  and  a  constant,  which  is  determined 
■Prom  the  experimental  data  ^on  the  confinement  casing  around  the 
charge10.  The  constant.,  A,  when  set  to  zero,  represents  an  unconfined 
explosive  charge  and  the  Defourneaux  relationship,  as  illustrated  in 
Figure  1.  <J>q  and  K  are  functions  of  the  angle  of  incidence,  i,  and 


9 E .  M.  Pugh ,  R.  J.  Eichelberger  and  N.  Rosvoker,  "Theory  of  Jet  Forma¬ 
tion  by  Charges  with  Lined  Conical  Cavities ,"  J.  Appl,  Physics ,  Vol. 
23,  No.  5,  May  1952. 

l0R.  DiPersio ,  J.  Simon ,  and  T.  Martin,  "A  Study  of  Jets  From  Scaled 
Conical  Shaped-Charge  Liners,"  BFLMR-1298,  August  1960.  (AD  #246332) 
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are  determined  for  certain  types  of  ewuosives.1.1  Figure  2  illus¬ 
trates  the  functional  relationship  between  4>0,  K,  and  i.  K  is  a 
constant  over  the  range  of  i  considered  for  typical  shaped-charge  de¬ 
signs,  i.e.  a  conical  liner  contained  in  a  cylinder  of  explosive. 
Figure  3  illustrates  the  linear  relationship,  Equation  1,  that  l/$ 
has  with  the  ratio  of  liner  mass  per  unit  volume  to  explosive  mass  per 
unit  volume,  given  as 


v  *■ 


p  c 

°H.E.e 
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where  pu  _  is  the  density  of  the  explosive.  The  values  of  l/$ 
n.u.  o 

and  K  used  in  Equation  1  are  the  Y-intercept  and  the  ratio  of  the 

slope  of  the  line  to  the  density  of  explosive  respectively.  Equation  1 

along  with  the  Taylor  formula, 

v  =  2D  sin  4  , 
o  l 

where  D  is  the  explosive  detonation  rate,  will  result  in  collapse 
velocities,  vq,  obtained  by  Gurney.12  Two  types  of  explosive  com¬ 
positions  are  shown  in  Figure  3  (data  taken  from  reference  11)  which 
represent  the  linear  function  at  a  constant  grazing  (parallel)  inci¬ 
dence,  i,  of  the  detonation  from  the  normal  to  the  metal  surface  (see 
Figure  1) . 

From  the  theory  of  Defoumeaux,  as  the  detonation  wave  sweeps 
toward  the  base  of  a  typical  shaped  charge,  $  decreases  due  to  the 
decrease  in  the  explosive  thickness,  e,  shown  in  Figure  3.  This 
assumption  that  <f>  decreased  monotonically  with  a  decrease  in  e  is 
justified  for  most  of  the  liner  collapse  since  there  is  sufficient 
time  for  the  liner  to  undergo  several  shock  reverberations  and 
achieve  a  bending  angle  close  to  its  maximum  before  entering  the  f^ow 
of  jet  formation.  However,  the  region  near  the  apex  of  the  cone 


1 1  /V.  DeFoumeaux  and  L.  Jacques ,  "Explosive  Deflection  of  a  Liner 
as  a  Diagnostic  of  Detonation  Flows,"  Proceedings  Fifth  Symposium 
(Interrtational)  on  Detonation,  ACR-184  Office  of  Naval  Research- 
Department  of  Navy ,  pp.  457-466 ,  Pasadena ,  California,  August  18~21, 
1970. 

12  1.  to.  Gurney ,  "The  Initial  Velocity  of  Fragments  from  Bombs,  Shells , 
and  Grenades,"  BRL  Report  No.  405,  Sept.  1943.  (AD  MTI36218) 


wiME.  ■C'MBi iisr 


TJEL. 


| 


i 

3 

4 


13 


t 


Figure  2 


0 


20  40 


60 


SO 


ANGLE  OF  INCLINATION  i  (degrees) 


4>o  and  K  are  functions  of  the  detonation  wave  angle 
to  the  liner. 
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Figure  3.  l/$  vs  y  as  determined  from  experiments.  Data 

is  taken  from  Reference  11. 


enters  the  flow  shortly  after  its  initial  acceleration  and  hence  does 
not  achieve  its  ultimate  bending  angle.  Therefore,  for  material  close 
to  the  liner’s  apex,  $  will  increase  to  a  maximum.  This  maximum  will 
be  located  near  liner  elements  that  originate  from  a  position  approxi¬ 
mately  40%  of  the  liner  height  measured  from  the  apex  of  the  cone. 

After  this  point,  (j>  will  decrerse  according  to  theory.  We  call  this 
"the  inverse  collapse  process."  Associated  with  this  process  is  the 
total  time  that  a  liner  element  takes  to  reach  the  axis.  The  equation 
for  the  collapse  time  is 

t  ■  fr  sin  i)/(D(sinCa+<j>)-sina)) ,  (2) 

where  several  new  variables  are  introduced,  r  is  the  instantaneous 
distance  of  a  liner  elemem  from  the  axis  before  its  collapse,  o 
is  the  half  angle  of  the  cone  of  a  conical  liner  or,  in  the  general 
case,  the  instantaneous  angle  with  the  axis  made  by  a  tangent  line 
to  the  liner,  t  will  be  affected  by  the  inverse  collapse  process. 
Material  close  to  the  liner  apex  enters  the  flow  of  jet  formation 
sooner  than  would  be  predicted  by  the  collapse  of  Defourneaux.4 

In  order  to  calculate  the  inverse  collapse  process,  a  new  equation 
for  the  bending  angle,  <j>  and  an  iterative  scheme  were  added  to  the 
code.  The  new  equation  for  4>  is 

4>n  a  4>  e  /~b_  ,  (3) 


where  b  =  Cj 


is  a  constant  determined  from  shape-charge  collapse  data.13  Th® 
iterative  scheme,  between  Equations  (2)  and  (3),  continues  until 
the  criteria  of  approaching  $  within  an  epsilon  is  satisfied. 

Having  determined  4>  and  t,  we  can  proceed  with  the  collapse  process. 


13F.  E.  Allison  and  R.  Vitali,  "An  Application  of  Jet  Formation  Theory 
to  the  105-rrtn  Shaped-Charge,"  BRLR-1165 ,  March  1962.  (AD  #277458) 
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The  velocity  of  collapse  of  the  liner  walls,  v  ,  toward  the 
charge  axis  to  font,-  the  jet  is  given  by 

v0-2D(sin  (*/2))/  fsin  i)  (4) 

The  apparent  explosive  detonation  velocity,  with  respect  to  liner 
wall  is  given  by 


Da»D/sin  i  (5) 

The  substitution  of  Equation  S  into  Equation  4  yields 

v0«2Da  Sin 0/2)  .  (6) 

This  is  the  so-called  Taylor  formula  utilized  in  the  code. 
The  angle  between  the  collapse  direction  and  the  charge  axis  is 
given  by 


y  -  0/2)  -  (a  +  0/2)).  (7) 

The  liner  element  first  hits  the  axis  at  a  distance,  sp,  from  the 
liner  apex,  given  by 

sp  =  (z  sin  <|>)  /  (cos  a(sin  (ot+<f>)  -  sina)),  (8) 

where  z  is  the  axial  component  of  the  liner  element  position  before 
its  collapse. 

While  the  liner  is  collapsing,  the  angle  formed  by  a  tangent 
to  the  collapsed  portion  on  the  axis  and  the  axis  itself  (called 
the  collapse  angle)  is  computed  from 


j 

i 


tan  fa  >  =  A2[sin  («+$)  -  sina]  tarnfr  +  rA<f>  cosot  (9) 
"  Azfsin  (a+<j>)~  -  sina]  -  rA<|>  cosot  tamf>  , 

where  Az  is  the  axial  increment  chosen  in  the  computational 
scheme,  and  A<J>  is  the  incremental  change  in  <j>  between  adjacent 
liner  elements.  The  cartesian  coordinates  of  a  liner  element 
(which  originates  at  position  z,r)  during  collapse  are  given  by 
the  pair  of  equations: 

x  (z,t)  =  z  +  v  nAt  sin  (a  +  (<J>/2)) 

0  (10) 
y  (z,t)  =  r  -  v^nAt  cos  (a  +  (4>/2)), 
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where  At  is  the  time  interval  taken  by  the  detonation  wave  between 
successive  liner  elements  and  n  is  a  positive  integer.  These 
equations  apply  in  the  time  interval  givon  by 

0  <  nAt  <_  t„ 

Figure  4  is  an  illustration  of  the  relationship  of  the 
variables  employed  in  the  BASC  code  of  a  generalized  axi symmetric 
collapse  of  a  typical  shaped  charge.  Shown  on  Figure  5  are  drawings 
giving  a  detailed  description  of  the  collapse  process.  Figure  5A 
shows  the  velocity  vectors  of  an  element  at  point,  P,  on  the  collapsing 
liner.  The  element  is  projected  toward  the  axis  of  symmetry  with  a 
collapse  velocity,  vq,  and  a  bending  angle,  When  the  detonation 

wave  with  velocity,  D,  has  progressed  a  distance,  P',  (i.e.  from  point 
P  to  point  Q)  during  the  time  interval,  t,  then  the  element  initally 
at  point  P  will  collide  with  the  cone  axis,  producing  the  geometrical 
relations  at  the  collision  or  stagnation  point,  spT  shown  in  Figure  5B. 
This  relationship  at  the  stagnation  point  is  with  respect  to  a  coordi¬ 
nate  system  moving  at  the  stagnation  point  velocity,  v_.  The  velocity 

sp 

of  the  liner  wall  flowing  into  the  stagnation  point  is  v^  and  the  angle 

between  it  and  the  cone  axis  is  the  collapse  angle,  3.  Figure  5C  shows 
a  cross-section  of  the  collapsing  liner  depicting  the  variables  employed. 
Applying  Bernoulli's  equations  at  the  stagnation  point,  we  find  that 
the  flow  velocity,  v^,  separates  into  two  equal  but  directionally 

opposite  velocities.  One  is  called  the  jet  velocity,  v.,  and  the  other 

J 

is  called  the  slug  velocity,  v^.  This  relationship  is  shown  on  Figure 

6.  Resolving  the  flow  velocity,  vf,  at  the  stagnation  point  in  the 

laboratory  coordinate  system  (Figure  6A) ,  the  following  set  of  equations 
are  obtained: 


v 


j 


+ 


v 


f 


(ID 

(12) 


In  accordance  with  the  laws  of  conservation  of  mass  and  momentum, 
when  the  liner  material  reaches  the  cone  axis  after  its  collapse,  it 
proceeds  either  as  a  fast-moving  jet  or  as  the  more  massive  but  slower- 
moving  slug  (Figure  6A) .  The  jet  velocity  equation  that  results  is 


Vj  »  Vq  cos (a  ♦  (*/2)  -  (6/2))/sin(S/2)  ,  (13) 

and  the  equation  for  the  slug  velocity  is 

vN  -  vq  sin(a  +  (*/2)  -  (3/2))/cos (3/2)  .  (14) 

The  relative  distribution  of  mass  (Figure  6B)  that  results  in  jet  and 
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A.  Velocity  vectors  of  an  element  of  the  collapsing  liner* 
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B.  Geometrical  relations  at  the  stagnation  point* 


Cono  Axis 


c.  Cross-section  of  collapsing  liner. 

Figure  5.  Illustration  of  the  relationship 
between  vaHables  and  the  collapse  process, 
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1  lustration  of  the  resolved  variables  in  the  laboratory 


slug  material  are  calculated,  respectively,  by 


***  2 

3m  “  sin 

and  (1 

■  cos2(3/2)  . 

Also,  the  relative  distribution  of  the  kinetic  energy,  for  the  jet 
and  slug,  respectively,  are 


3^“  COS2 (a  +  U/2)  -  (3/2)) 


ocN  2 

a  sin  (a  +  C4>/2) 


(3/2))  . 


The  variables  dm^  and  dE^  are  the  incremental  change  in  the  liner's 
mass  and  kinetic  energy,  respectively. 

The  impinging  flow  velocity,  v  fl  hat  been  considered  by  shaped- 

^  14  15 

charge  researchers  such  as  Walsh,  ev\  al,  ;  Cowan,  et.  al.  ;  and 
Chou,  et.  al.8  to  be  critical  in  the  jet-formation  theory.  As  illus¬ 
trated  in  Figure  7,  when  vf  is  less  than  the  material  sound  speed,  c, 

the  jet  formed  is  coherent  or  a  good  jet  (Figure  7A) .  Even  whenvf 

is  slightly  greater  than  c,  this  too  forms  a  coherent  jet  (Figure  7B) . 
But,  when  vf  is  sufficiently  greater  than  c,  the  jet  will  be  incoherent 

or  bifurcated.  Wc  call  it  a  no-jet  condition  (Figure  7C).  From 
equations  11  and  12,  solving  for  v^  yields 

vf  a  -5(vj  -  vtP  •  (17) 

We  then  use  the  following  relationship  as  the  jet  limiting  criteria 
for  a  coherent  jet: 

Mr  s  vf/<i  <  1,23  •  (IB) 


V.M.  Walsh,  R.G.  Shreffler ,  and  F.J.  Willig ,  "Limiting  Velocity  Con¬ 
ditions  for  Jet  Formation  tn  High  Velocity  Collisions ,"  Journal  of 
Applied  Physics,  Vol.  24,  No.  3,  pp.  349-359,  March  1957. 

ISG.R.  Cowan  and  A.H.  Holtzman,  "Flow  Configurations  in  Colliding 
Plates:  Explosive  Bonding ,"  Journal  of  Applied  Physics,  Vol.  34, 

No.  4,  pp.  928-939,  April  1963. 
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Figure  7.  Illustration  of  the  variables  employed  in  the  BASC  code 
to  limit  jet  formation.  The  jet  limiting  criteria  for 
a  coherent  jet  (A  and  B)  is  MR  =  v^/c  <  1.23. 


The  threshold  factor  for  a  coherent  jet,  KL-1.23,  is  determined  by 
comparing  rtuiographic,  jet -particle  data  measured16  from  experi¬ 
ments  with  BASC  results.  The  threshold  factor  holds  for  several 
materials  considered  in  BASC  calculations. 


In  order  to  determine  the  jet-tip  velocity  and  the  mass  of  the 
lead  pellet  to  be  used  in  the  shaped-charge  penetration  theory,  the 
inverse  velocity  gradient  has  to  be  equilibrated  and  the  jet  mass  in 
this  region  to  be  compresses  into  a  single  zonal  element,  LP.  This 
zone  will  then  contain  the  so-called  steady-state,  lead  pellet.  The 
equilibrated  jet-tip  velocity,  v.  ,  is  obtained  by  the  following 
process: 


,  i+l  1  th 

for  Vj  >  Vj  and  1  <_  i  <_  LP,  where  i  is  the  i  zonal  element.  The 


equilibrated  jet-tip  velocity  is 


(18) 


The  steady-state. 


lead  pellet  mass  is 


ii 


(19) 


In  the  theory  of  shaped-charge  jet  penetration  into  a  target 
used  in  BASC,  an  important  parameter  is  the  time  that  a  given  liner 
element,  which  enters  into  the  jet,  arrives  at  the  bottom  of  the 
target  hole  when  penetration  is  in  progress.  Time  is  usually  started 
from  the  moment  the  detonation  wave  reaches  the  apex  of  the  liner. 

A  time  parameter,  9,  is  defined  by 


e  *  tz  +  T  ,  (20) 


16  Private  Cormuniaation  from  J.  Blisohe  at  BRL. 

l7R.  DiPcrsio ,  W.H,  Jones ,  A.B.  Merendino ,  and  J.  Simon ,  " Character - 
istio8  of  Jets  from  Small  Caliber  Shaped  Charges  with  Copper 
and  Aluminum  Liners,"  BRLMR  No.  1866 ,  September  1967.  (AD  # 823839 ) 
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where  tz  is  the  time  at  which  the  detonation  wave  reaches  the  liner 

element  at  an  initial  axial  distance,  z,  from  the  apex  and  t  is  the 
time  taken  by  this  element  to  collapse  onto  the  charge  axis.  The 
sum  of  standoff  distance  from  the  liner  base  to  the  initial  surface 
of  the  target,  h,  and  the  liner  height,  H,  defines  Zq: 

lQ  *  h+H.  (21) 

The  time  for  the  tip  of  the  jet  to  reach  the  target  surface  is 

To  ■  VvJo  C22) 

It  is  assumed  that  the  jet  tip  originates  from  the  zonal  element,  LP, 
and  contains  the  highest  velocity  of  the  jet,  Vj0,  An  element  from 
the  liner  between  LP  and  base  results  in  a  jet  velocity,  v^ ,  which  is 
less  than  v.  and  is  a  function  of  its  initial  axial  distance,  z.  A 
parameter,  y  is  defined  by 


u  =  v.0/Vj  US) 

The  time  for  the  portion  of  jet  formed  from  a  given  element  of  the 
liner  to  reach  the  bottom  of  the  target  hole  is 


T  =  y(1+k)  [T  +  (1+k)  |  y~k (Af ) . ]  (24) 

O  0  1 

where  (Af) .  *  f.  -  f.  ,  and  denotes  the  i  n  zonal  element.  In 

this  equation,  k  =  /p j7pc »  where  p.  is  the  jet  density  (assumed 
equal  to  the  liner  density),  and  Pc Jis  the  target  density,  f  is  a 
time  parameter  given  by 

f  =  ((0/p)  -  (ip/vj0))  ,  (25) 

where  "sp  was  previously  defined  by  Equation  8,  and  v.  ,0,  and  y 
were  defined  by  Equations  18,  20,  and  23.  The  value^8f  u  to  be 
used  outside  the  bracket  in  Equation  24  is  that  which  applies  to 
the  original  elemental  liner  position.  Af  is  the  difference  in  f 
values  between  adjacent  elemental  positions,  and  the  summation 
applies  to  all  adjacent  elemental  positions  up  to  the  point  in 
question  on  the  liner. 


While  the  jet  is  still  continuous  in  nature,  it  stretches  due 
to  its  velocity  gradient  and,  therefore,  decreases  in  diameter  with 
increase  in  time.  The  equation  for  calculating  the  jet  radius,  r. 
is 


r 
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2  re  D  sin 

TT 
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fir51" 


In  this  equation,  r  is  the  initial  radial  position  of  the  liner 
element  from  the  charge  axis,  e  is  the  liner  thickness  at  this 
position,  and  6  is  the  collapse  angle  that  is  formed  when  this 

element  reaches  the  charge  axis.  The  factor, ,  is  given  by 

0  z  T 


AZ 

At 


z  T 


AVj 

At 


v  +  A(z-r  cot(B/2)) 


(27) 


In  this  equation,  t  is  defined  following  Equation  20  and  At  is 
the  incremental  time  interval  between  arrival  of  the  detonation 
wave  between  successive  liner  elements  (a  constant) .  The  value 
computed  by  Equation  27  is  a  function  of  time,  T,  which  starts 
at  zero  when  the  detonation  wave  first  reaches  the  liner  apex.  It 
is  a  negative  value  which  increases  in  absolute  value  as  T  increases. 
Therefore,  it  can  be  seen  from  Equation  26  that  the  jet  radius,  which 
originates  from  material  at  any  point  on  the  liner,  decreases  with  an 
increase  in  T.  The  minus  sign  in  Equation  26  is  necessary  to  make 
r.  a  positive  quantity. 

When  the  jet  cannot  sustain  any  further  stretching,  it  breaks 
up  into  individual  axial  particles.  It  is  assumed  that  this  occurs 
throughout  the  whole  jet  at  the  same  time.  The  breakup  time  fox 
the  jet  is  designated  as  T.  and,  at  present,  must  be  obtained  from 
experimental  observations.  At  times  greater  than  the  jet  breakup 
time,  the  individual  jet  particles  do  not  stretch  in  length  or 
decrease  in  diameter.  The  constant  jet  particle  radius  that  one 
obtains  for  material  originating  from  a  given  element  of  the  liner 

AZ 

is  calculated  by  Equation  26  in  which  the  factor  Tr—|>r  is  calculated 

z 

at  the  breakup  time,  T  in  Equation  27.  However,  radii  of  different 
particles  are  different  due  to  the  variability  of  the  initial  liner 
element  position  in  Equations  26  and  27. 

The  equations  that  are  used  for  jet  penetration  theory  are 
dependent  upon  the  standoff  distance  between  the  charge  and  the 
target.  If  the  target  is  placed  close  enough  to  the  charge  so 
that  the  jet  tip  reaches  the  target  before  the  time  of  jet  breakup, 
one  set  of  equations  are  used.  On  the  other  hand,  if  the  jet 
particulates  before  reaching  the  target  surface,  a  different  set 
of  equations  must  be  used.  In  the  former  case,  even  though  the  jet 
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starts  penetrating  in  the  continuous  state,  it  becomes  discontinuous 
before  the  end  of  penetration  and  different  equations  are  required 
after  time  T^.  The  variables  and  the  two  states  of  penetration  are 
illustrated  in  Figure  8. 

The  jet  penetrates  the  target  in  both  the  continuous  and  discon 
tinuous  state  whenever  the  following  criteria  is  satisfied: 

Z  <  v.  T,  .  (28) 

o  jo  1 


The  depth  of  penetration  into  the  target,  p,  for  any  T,  such  that 

T  <  T  <  T. , 
o  —  —  1 


I>  .  2  [(T_,  (*/(*.») 

0 


(29) 


where  all  factors  have  been  previously  defined.  For  times  greater 
than  T,  the  equation  used  is 


p  =  ZQ  [(k+1)  ,  T 


TA  (k/ (k+1))  T 


Tj+kTj  “ 


(30) 


The  total  penetration  depth  into  the  target  is  calculated  by 


k 


PT  ’ 


V 


k_  _J_ 

k(k+l)UminTt (v .  T.)\+1Z  k+1-Z 
v  '  r  jo  1  o  o 


(31) 


The  only  hitherto  undefined  term  in  Equation  31  is  this  factor  U 
This  is  called  the  minimum  penetration  velocity.  Accoiding  to 
theory,  the  velocity  of  penetration  into  the  target  mono  Ironically 
decreases  with  increase  in  penetration  depth  until  it  reaches  the 

value  Umin.  When  Umin  is  reached,  penetration  stops  and  all' 
remaining  jet  material  just  piles  up  at  the  bottom  of  the  target 

hole.  Umin  is  an  empirical  constant  whose  value  depends  upon  the\ 
target  material  and  its  hardness.  It  is  invariant  with  standoff 
distance.  With  a  given  charge,  target,  and  standoff  distance,  the 
computer  calculates  the  constant  total  penetration  depth  from 
Equation  31.  This  value  is  used  by  the  computer  as  a  signal  to 
stop  calculating  p  in  Equation  30  and  also  to  determine  the  time 
at  the  end  of  the  penetration  process.  The  radius  of  the  hole  in 
the  target  before  the  time  of  jet  breakup  is  given  by 
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Figure  8. 


:hematic  cross  section  of  »  sh.pcd-ch.rged  jet 
.netratine  a  target  in  two  states. 


(a)  Continuous  and  discontinuous, 


(b)  Discontinuous. 
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This  equation  contains  one  J«st  empirical  constant  of  the  target, 
ck*  This  is  the  jet  kinetic  energy  needed  to  produce  a  unit  of  hole 


volume  in  the  target.  Its  value  depends  upon  the  target  used  and 

?fiSailHr77eSSl^The  ualv®  °f  ri  in  E£luation  32  is  obtained  by  Equations 
26  and  27.  After  the  breakup^™.  T],  the  hole  radius  is  obtained  by 
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v.  T, 

JO  1 
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(33) 


where  r.  is  now  obtained  by  Equations  26  and  27  with  the  restriction 
AZ 


that  the  factor 


At 


It  must  be  evaluated  at  T  =  T^.  The  hole  profile. 


as  computed  by  Equation  33,  is  terminated  when  the  penetration  depth, 
p,  reaches  the  value  P^,  as  calculated  by  Equation  31. 


The  last  set  of  penetration  equations  is  used  when  the  standoff 
distance  between  charge  and  target  is  so  large  that  the  following 
condition  is  satisfied: 


Zo  >  Vjo  T1  ‘  (34) 

In  this  case,  all  penetration  is  accomplished  by  the  jet  while  it 
is  particulated  in  nature.  The  equation  for  penetration  depth  is 
then  r 


p  =  v.Q  (T-To)  Tj  /  (Tj  +  CT/k))  (35) 

where  the  time  factor,  T, “Varies  between  the  time  of  first  target 
impact,  To,  and  the  time  of  last  jet  penetration,  T  .  The  total 

penetration  depth  is  given  by 

't  ■  fvjo  T1  -  /il""  1,  (vjo  Tj  *  (2o/k)]k 


(36) 
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The  radius  of  the  hole  produced  in  the  target  as  a  function  of  its 
depth  is  given  by 

rc  -  vjQ  rj  Cl-Cp/k  vjQ  T^)  /  p^/ik  c^  ,  (37) 

where  the  jet  radius,  r j ,  is  obtained  by  Equations  26  and  27  with 

the  restriction  that  the  factor  bl  i  is  evaluated  at  T  *  T.. 

At  >T 
z 

The  variation  of  total  penetration  depth  with  standoff  distance 
is  computed  by  means  of  Equations  31  and  36,  The  factor  Zq  is  the 

only  variable  in  these  equations.  Equation  31  is  used  first  until 
Z  increases  to  the  value  v.  T, ,  then  Equation  36  is  used. 

O  * 

III.  CALCULATIONAL  SCHEME 

The  BASC  code  enables  one  to  perform  parametric  studies  for 
designing  warheads.  The  variables  employed  in  the  code  of  the 
generalized,  axisymmetric  collapse  of  a  shaped  charge  were  illustrated 
previously  in  Figure  4.  The  parameters  which  can  vary  include  the 
following: 


a 

K 

E 

P. 


H 


(ALPHA) 

The  half  angle  of  the  liner  (degrees) 

(CON) 

The  empirical  constant  for  the  detonation 
products.  Value  known  for  Composition  B 
explosive. 

(EPS) 

The  thickness  of  the  liner 

(cm) 

(RHOJ) 

The  density  of  the  liner 

(gm/cm3) 

(RHOC) 

The  density  of  the  targe 

(gm/cm3) 

(RF) 

Radius  of  the  base  of 

targe  (cm) . 

(H) 

Height  of  the  liner 

•f  H  is  zero, 

H  will  be  calculated  by  H 

rF  . 
tana 

30 

D 

(D) 

Detonation  velocity  (cm/ys) 

SO 

(SO) 

Standoff  distance  between  base  of  liner 
and  target,  (cm) 

CK 

(CK) 

Constant  for  determining  hole  volume 

T1 

(Tl) 

Jet  breakup  time  (ys) 

if™ 

(UMIN) 

Velocity  cutoff  for  the  penetration  of  jet 
into  a  target  (cm/ys) 

L 

(DPOINT) 

Initially  the  total  length  of  the  charge 
(cm)  but  is  converted  to  be  the  initiation 
point  of  the  explosive  (i.e.  DPOINT  =  L-H) 

RDPT 

(RDPT) 

Radius  above  axis  where  explosive  is 
initiated  (cm) . 

JOHN 

(JOHN I) 

If  JOHNI  greater  than  zero,  <J>  will  vary; 

if  not.,  $  will  be  a  calculated  constant. 

0 

N 

(N) 

The  number  of  zones  into  which  the  grid 
is  subdivided.  If  N  is  zero,  default  is 
seventy  zones. 

The  code  will  set  up  the  grid  based  upon  FiguT’s  4  and  the 
equation  for  increment  in  Z  is 

DZ  =  H/N  (38) 

and  from  this  the  time  increment  is 

DTZ  =  DZ/D  .  (39) 

Therefore,  at  each  increment  in  Z  we  have  also  the  corresponding 
time  increment. 


The  code  then  marches  sequentually  through  the  governing  equa¬ 
tions  (1-19),  calculating  and  storing  jet  formation  information  to 
be  printed  and  plotted  at  the  end  of  the  iteration.  Those  vari¬ 
ables  include:  , 

•  r  ,th 
i  I  i  zone 


ct 


z 


Al.PHA (I)  Tangent  angle  of  liner  wall  to  the  axis 
(degrees) 


Z(I)  z  distance  from  apex  to  base  of  liner 
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1 


I 

.  t  t 

<  Z 

TZ(I) 

tih 

Time  at  the  i  zone 

1  , 

l  Y 

GAMMA(I) 

Angle  (degrees) 

1 

J 

i 

i 

E 

E(I) 

Thickness  of  explosive  above  i**1  zone 

l 

1 

\ 

$ 

PHI  CD 

Angle  (degrees) 

1 

e 

BETA (I) 

Angle  (degrees) 

1 

A4> 

DPHI(I) 

Increment  of  angle  <f> 

!  1 

I  ♦ 

RPHI(I) 

Reciprocal  of  angle  4> 

! 

V(I) 

Collapse  velocity  (cm/ys) 

■  ]  r 

R(I) 

Radius  of  the  i^  zone  (cm) 

1 

•  T 

TAU(I) 

Time  to  collapse  i^  zone  (ys) 

j  sp 

CCD 

Distance  from  the  apex  of  the  liner  where 

ith  zone  will  hit  the  axis  (cm) 

vjcn 

Velocity  of  the  jet  in  (cm/us) 

1  V 

;  VN 

« 

VN(I) 

Velocity  of  the  slug  in  (cm/ys) 

dm. 

.  ast 

DMJ(I) 

Relative  inss  of  the  jet,  dimensionless 

dy 

dmL 

DMN(I) 

Relative  mass  of  the  slug,  dimensionless 

dE. 

DEJ(I) 

Relative  energy  of  the  jet,  dimensionless 

i 

j  . 

;  dEL 

dEN 

DEN  ( I ) 

Relative  energy  of  the  slug,  dimensionless 

r  • 

dEL 

P  "'  N 


‘I 


,1 

fc  - 

i 

4 


RV  Cl) 


Relative  velocity  between  the  jet  and  slug 
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In  the  penetration  portion  of  BASC,  the  constants  listed  below 
must  be  calculated: 

VJO  =  jet  tip  velocity  from  Equation  18  (40) 

ZO  =  H  +  SO  (41) 


PT  =  Total  penetration  (cm)  from  Equation 

31  or  36  (44) 

The  next  set  of  outputs  are  for  the  penetration  of  the  jet 
into  the  target.  The  code  also  marches  through  the  governing 
Equations  20-37,  calculating  and  storing  penetration  information 
to  be  printed  and  plotted  at  the  end  of  the  iteration.  The  vari¬ 
ables  calculated  and  stored  are  listed  below: 


i 

I 

.th 

i  zone 

U 

AMU(I) 

Relative  velocity  between  the  jet  tip 
velocity  and  jet  velocity  of  the  i  zone 

G 

THETA(I) 

Time  parameter  in  microseconds 

f 

F(I) 

Time  parameter  in  microseconds 

Af 

DF(I) 

Time  increment  of  f  in  microseconds. 

T 

TCI) 

t.h 

Time  that  the  i  element  reaches  the  bottom 
of  the  target  hole  in  microseconds. 

At 

DT(I) 

Time  increment  of  T 

t 

RSQ(I) 

is  ^r j 2  which  is  the  radius  of  the  it*1 
element  of  the  jet. 

A 

A(I) 

is  equal  to  z  -  r  cot  p/2. 

AA 

Av . 

DELA(I) 

DVJ(I) 

is  the  increment  of  A(I). 

is  the  increment  of  jet  velocity. 
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DZODT(I)  Equation  27  in  governing  equations 


rc  RC(I)  Radius  of  the  hole  in  target  of  the  ith  zoi 

in  centimeters 

p  p(I)  Depth  of  penetration  in  centimeters 

After  these  parameters  have  been  printed  and  plotted,  the  code  then 
returns  to  start  and  begins  another  case.  This  is  continued  until 
the  end  of  file  (i.e.  next  problem)  is  encountered  causing  the 
run  to  terminate. 


IV.  COMPARISON  OP  BASC  CODE  RESULTS 

The  performance  of  the  BASC  code  is  best  illustrated  by  pre¬ 
senting  results  of  a  calculation  and  comparing  these  results,  where 
possible,  with  experimental  data  or  with  results  from  oilier  numeri¬ 
cal  techniques.  The  first  set  of  comparisons  will  be  with  experi¬ 
mental  results  from  the  following: 

a.  105-mm,  unconfined,  4'*°,  copper-lined  shaped  charge 
with  a  Composition  B  explosive  fill  tested  by  Allison 
and  Vital!13. 

b.  A  study  of  jets  from  scaled,  42°,  copper-lined,  conical 
shaped  charges  filled  with  Composition  B  explosive  (test 
by  DiPorsio,  et.  al.10). 

c.  3,g-inch,  BRL  precision,  shaped  charge  with  a  copper  liner 
42°  apex  angle,  and  Composition  B  explosive  fill,  having  its 
jet  characteristics  measured16  from  radiographic  data 
recorded  at  the  BRL. 

d.  A  study  of  jet  characteristics  from  small-caliber  shaped 
charges  with  copper  and  aluminum  liners  and  variation  in 
apex  angle  from  20°  to  90°.  All  charges  were  filled  with 
Composition  B  explosive.  The  tests  conducted  at  the  BRL 
by  DiPersio,  et.  al?7 

The  second  set  of  comparisons  will  be  another  numerical  technique. 

The  other  technique  is  the  two-dimensional,  hydrodynamic  computer 
program  based  upon  the  Eulerian  numerical  scheme  called  the  BRLSC 
(Ballistic  Research  Laboratory  Shaped  Charge)  code.18  The  BRLSC 

16W.L.  Gittinga,  "BRLSC:  An  Advanced  Eulerian  Code  for  Predicting 
Shaped  Charges >"  Vol ,  I,  BRL  CR  279 ,  Prepared  by  Syatemt  Science 
and  Software,  Deaember  1978.  (AD  M02S9G2) 
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code  is  a  modified  version  of  the  HELP  (Hydrodynamic  Elastic-Plastic-) 
code19  developed  by  System,  Science  and  Software. 

Figure  9  and  10  are  experimental  comparisons  for  the  105-mm, 
unconfined,  shaped  charge.  Collapse  velocity,  v  ,  and  jet  velocity, 
vj  as  a  function  of  the  relative  distance  from  the  apex  of  the  cone 
are  shown  on  Figure  9.  The  dashed  line  is  the  jet  velocity  after  jet 
compression  and  illustrates  exact  agreement  with  both  the  jet  tip  velo¬ 
city  and  position  on  the  axis  between  BASC  and  experiment.  This  illus¬ 
trates  that  jet  tip  is  compressed  into  one  element  at  a  position  which 
is  approximately  40%  of  the  distance  from  the  apex  of  the  cone.  The 
collapse  velocity  is  slightly  higher  than  that  calculated  by  Allison 
and  Vitali,  but  the  overall  agreement  is  good.  Figure  10  shows  a 
cjmparison  of  the  predictions  of  jet  mass  versus  cone  mass  between 
the  BASC  code  and  Allison  and  Vitali.  Those  two  predictions  are  in 
slight  disagreement  with  one  another  because  Allison  and  Vitali  were 
unable  to  recover  100%  of  all  the  slug  material.  They  used  essentially 
the  same  theory  as  the  BASC  code  to  predict  jet  mass,  and  all  of  the 
material  is  required  to  accurately  predict  the  true  conservation  laws. 

Figure  11  is  the  second  experimental  comparison  of  the  scaled, 
shaped  charge  with  a  heavy  confinement  casing.  This  figure  shows 
collapse  velocity,  v  ,  and  jet  velocity,  v- ,  as  a  function  of  the 
relative  distance  fr8m  the  apex  of  the  ccne.  The  dashed  line  is  the 
compressed  jet  velocity  distribution.  This  again  illustrates  exact 
agreement  for  both  jet  tip  velocity  and  its  positior  on  the  axis 
between  BASC  and  exrdriment.  The  open  squares  nt  approximately  the 
48%  position  fr~’M  the  apex  of  the  cone  shows  the  spread  in  the  experi¬ 
mental  data  from  the  experimental  scaled  rounds.  The  BASC  code  results 
arc  identical  for  the  same  scaled  rounds  (see  Reference  10).  Shown 
at  the  base  of  the  cone,  for  values  greater  than  80%  of  the  distance 
from  the  apex  of  the  cone,  is  a  change  In  the  slope  of  the  jet  velo¬ 
city  curve.  This  is  due  to  gas  leakage  or  breaking  of  the  confinement 
casing.  This  phenomenon  is  modeled  in  BASC  as  a  gradual  change  in  the 
confinement  factor,  A  in  Equation  1,  until  it  reaches  zero,  i.e. 
unconfined.  This  comparison  shows  very  good  agreement  with  the  experi¬ 
mental  results. 

Table  I  is  a  tabulation  of  a  comparison  between  measured, 
radiographic  experimental  data  and  BASC  code  results  from  the  3.2- 
inch,  BRL,  precision  shaped  charge.  The  jet  tip  velocity  and  mass 
at  jet  particle  number  one  are  in  agreement,  but  the  accumulated 
total  jet  mass  from  jet  particle  number  one  to  jet  particle  number 


101 Hageman  and  J.M.  Walah ,  " HELP ,  A  Multi-Material  Eulericm 
Program  For  Compressible  Fluid  and  Elastic-Plastic  Flows  in  Two 
Space  Dimensions  o'ld  Timet"  BRL  CR  39,  Prepared  by  System ,  Science 
and  Software ,  May  1971.  (AD  Nos.  726459  and  726460) 
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Figure  9.  Comparison  between  experimental  and  BASC  Code  results 
of  jet  and  collapse  velocities  from  the  105-mm,  ancon 
fined,  shaped  chavge. 


36 


ML  in  /£•  ii,t it  A.  i-  ;  / 


L  l‘-1  a  ^ 


± 

« 

M 


!-■  ,  I _ I _ L 

3  £  2 

6 *ssvw  lar  - 


o 

♦ 


37 


VJ  (MM/MICROSEC) 


- 


(l-H)  =  20. 5537  MM 

RD  s  0.0000  MM 
EXPERIMENTAL  DftTfi  «. 


0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8'  0.9 

RELATIVE  DISTRNCE  FROM  CONE  APEX 

filPHR  s  21.00  DEB 
RF  =  24.0030  MM 

EPS  =  1.78459  MM 


Figure  11.  Comparison  between  BASC  coda  and  experimental  results 

for  scaled,  heavily-confined,  shaped  charges  (Reference 
10).  Jet  and  collapse  velocities  vs  %  of  the  distance 
from  the  apex  of  the  cones  are  shown.  Solid  lines  are 
BASC  code  results.  The  dashed  line  is  the  final,  computed, 
compressed  jet  velocity  profile.  The  jet  tip  for  both 
experimental  and  BASC  results  is  shown  compressed  into 
a  position  48%  from  the  apex  of  the  cone. 


forty- three  is  10  grams  heavier  in  the  BASC  code  results.  Also,  the 
code  results  show  46%  more  total  accumulated  kinetic  energy  from  the 
first  forty- three  jet  particles  than  that  obtained  from  experiment. 
Overall,  this  comparison  is  reasonable  considering  the  material  loss 
from  such  physical  effects  as  erosion  and  possible  errors  in  the  data 
reduction  measurements  and  calculations  of  the  experimental  results. 


Table  I.  BASC  Calculations  and  Measured,  Experimental  Radio- 
graphic  Data  of  Jet  Velocity,  Mass,  and  Kinetic 
Energy  for  the  3.2-in.  (81.3-mm),  BRL,  Precision, 

Shaped  Charge 

jet  Velocity  (km/s)  Accumulated  Mass  (g)  Kinetic  Energy  (kJ) 


Particle  No. 

Measured 

BASC 

Measured 

BASC 

Measured 

BASC 

1 

7.74 

7.74 

3.7 

4.3 

110. 

128. 

43 

2.90 

2.98 

22.6 

32.6 

348. 

507. 

The  last  experimental  comparison  will  be  the  characteristics  of 
jets  from  small  caliber  shaped  charges  with  copper  and  aluminum  liners. 

In  this  study,  only  the  liner  designs  referred  to  as  li-inch  (38.1-mm) 
liner  base  diameter  chargeg  in  Reference  17  will  be  considered.  The 
apex  angle  will  include  20°,  40°,  60°  and  90°  cones  for  both  the  cop¬ 
per  and  aluminum  liners.  These  are  all  confined  with  a  steel  casing 
around  the  charge.  The  pertinent  dimensions  for  the  shaped-charge 
designs  considered  can  be  found  in  Reference  17  on  pages  9  and  10.  The 
results  from  the  BASC  code  calculations  and  the  experiments  are  sum¬ 
marized  on  Table  II.  All  of  the  results  of  the  BASC  code  are  within 
a  5%  error  boundary  except  the  liner  type  with  a  90°apex  angle  and 
aluminum  cone  which  is  10%.  The  results  of  the  calculations  are  slightly 
higher  than  the  experiments  for  all  of  the  aluminum  liner  types. 


Table  II.  Summary  of  Results  From  the  Small-Caliber,  Shaped- 

Charge  Study  for  the  38.1-mm  (li-in.),  Base  Diameter 
Design 

Jet  Tip 

Velocity,  V^o  (km/s) 

Liner  Type  Experiment  Predicted  (BASC)  \ 


o 

o 

CM 

apex  angle,  Cu 

Cone 

9.9 

9.80 

O 

O 

apex  angle,  Cu 

Cone 

8.2 

8.00 

O' 

o 

o 

apex  angle,  Cu 

Cone 

6.7 

6.64 

90° 

apex  angle,  Cu 

Corie 

5.5 

5.48 

39 


JC.i/iiA.Wa  a!. 45:..  t. 


•  •  * 


Experiment 


20  apex  angle,  A1  Cone 
40°  apex  angle,  A1  Cone 
60°  apex  angle,  A1  Cone 
90°  apex  angle,  A1  Cone 


Predicted  (BASC) 

11.80 

9.75 

8.54 

7.50 


The  final  two  sets  ot  comparisons  involve  the  two  numerical  tech¬ 
niques:  BASC  and  BRLSC.  First,  Figure  12,  is  a  comparison  of  the 
calculated  flow  field  for  the  105-mm,  unconfined  shaped  charge  reported2 
for  the  BRLSC  code  and  the  BASC  code  at  25  ysec  after  initiation  of  the 
explosive.  The  radius  of  the  jet  and  slug  are  in  excellant  agreement, 
but  the  tip  of  the  jet  is  slightly  behind  in  the  BASC  calculation. 

This  is  due  to  the  fact  that  the  BASC  code  takes  very  large  time  steps 
which  will  not  allow  the  inverse  velocity  gradient  to  be  equilibrated 
at  this  snap  shot  in  time.  The  second  set  of  comparisons  between  BASC 
and  BRLSC  codes  is  shown  on  Table  III.  The  results  shown  are  from  an 
improved  version  of  the  BRLSC  Code,2^  Table  III  is  a  tabulation  of 
six  calculated  results  from  BASC  and  BRLSC  as  well  as  four  experimental 
results,  which  can  be  used  as  a  guide  in  this  comparison.  All  of  the 
calculations  and  experiments  in  this  study  involved  identical  copper 
liners  (21  ,  81.3-mm  base  diameter,  1.9-mm  thickness).  Calculations 
1  through  4  omployed  TNT,  Comp  B,  Octol  and  PBX,  respectively  as  ex¬ 
plosive  filler. 


Table  III.  Computational  Matrix  and  Summary  of  Results 


Calculation 
Numb  er 


HE 

Confinement 

Total  Jet 

Mass  (g) 

Jet  Kinetic 

Enerj 

'Fill 

BRLSC 

BASC 

BRLSC 

BASC 

TNT 

Light 

19.66 

28.42 

196. 

385. 

Comp  B 

Light 

22.56 

30.73 

287. 

507. 

Octol 

Light 

27.40 

32.03 

349. 

570. 

PBX 

Light 

27.97 

33.49 

381. 

614. 

Comp  B 

Heavy 

33.63 

60.41 

475. 

1262. 

Octol 

Heavy 

35.17 

60.16 

360. 

1406. 

20i?.  T.  Sedgwick ,  L.J.  Walsh  and  M.S,  Chawla,  "Effects  of  High  Explosive 
Parame bens  and  Degree  of  Confinement  on  Jets  from  Lined  Shaped 
Charges , "  BEL  CE  P.45,  Prepared  by  System ,  Saienoe  and  Software , 

July  1975.  (AD  ttB00698?L) 
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Calculation 
,  Number 

Predicted  Final 

Jet  Tip  Velocity  (km/s) 

Measured 

Jet  Tip  Velocity (km/s) 

BRLSC 

BASC 

1 

5.95 

6.79 

6.8 

2 

7.15 

7.74 

•  7.7 

3 

7.40 

8.17 

8.1 

4 

7.80 

8.44 

8.3 

5 

7.44 

8.86 

— 

6 

7.78 

9.36 

— 

In  these  four  calculations,  the  charge  was  confined  in  an 
aluminum  casing  (treated  as  being  unconfined  in  the  BASC  code) . 
Calculations  5  and  6  employed  Comp  B  and  Octol  as  explosive  fills 
which  were  confined  in  a  steel  casin  \  with  a  thickness  of  10mm. 

The  calculational  matrix  is  shown  in  Table  III  with  the  degree  of 
confinement  provided  by  the  aluminum  and  steel  casing  is  referred 
to  as  light  and  heavy,  respectively.  The  summary  of  the  results 
of  both  codes  (BASC  and  BRLSC)  as  well  as  the  results  of  the  jet 
Lip  velocities  measured  from  experiments21  are  also  shown  in  Tuble 
III.  The  quantities  summarized  represent  the  total  jet  which  is 
composed  of  jet  material  having  a  velocity  greater  than  or  equal  to 
a  velocity  of  2.8  km/s.  The  results  indicate  the  ratio  of  predicted 
total  jet  mass  of  BRLSC  code  to  the  BASC  code  is  approximately  65% 
for  all  calculations. 

The  comparison  of  jet  tip  velocity  data  is  illustrated  on 
Figure  13.  The  theoretically  predicted  values  for  calculations  1 
through  4  from  both  the  BASC  and  BRLSC  codes  are  compared  with  the 
experimental  determined  jet  tip  velocity  data.  Figure  13  is  a 
plot  of  calculated  jet  velocites  versus  the  measured  jet  tip  velo¬ 
cities.  The  open  triangles  are  the  final  jet  tip  velocities  as 
predicted  by  the  BRLSC  code20  and  the  open  squares  are  those  predicted 
by  the  BASC  code.  There  is  very  good  agreement  between  the  BASC  code 
and  the  observed  results  for  those  rounds  considered  and  the  BRLSC 
code's  agreement  is  less  than  10%  except  for  TNT  which  is  12.5%. 

In  summary,  we  have  demonstrated  in  these  comparisons  the  predic¬ 
tive  ability  of  the  BASC  code.  In  addition,  several  of  its  salient 
features  has  been  shown.  These  features  include  the  predictive  char¬ 
acteristics  of  BASC  with  respect  to  variations  in  casing  confinement, 
variations  in  cone  apex  angle,  variations  in  liner  density  and  thick¬ 
ness,  variations  in  explosive  fill,  and  the  effects  of  scaling.  We 
have,  also,  demonstrated  its  predictive  ability  with  respect  to  other 
sophisticated  numerical  schemes. 


2 V.  i Simon,  "The  Effect  of  Explosive  Detonation  Characteristics  on 
Shaped-Charge  Performance BRLMR-2414,  August  1979.  (AD  #B000337L) 
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Figure  12,  Comparison  of  calculations  of  the  flow  field  between 
BASC  code  ard  a  2-D,  Hydrodynamic  Code  (BRLSC") 
(Reference  2) .  These  calculations  are  of  the  105-mm, 
unconfined,  shaped-charge  shown  at  the  time  approxi¬ 
mately  25ps  after  the  initiation  of  the  explosive. 

The  BRLSC  Code  results  are  to  the  left  of  the  axis 
of  symmetry  and  the  BASC  code  results  are  to  the 

Tight , 


Based  upon  these  comp'risons,  it  is  concluded  that,  for  many  'oni- 
cal  shaped-charge  designs,  the  BASC  code  is,  a  useful,  predictive 
model  for  jet  characteristics.  Since  it  requires  only  fractions  of  a 
minute  tc  do  several  shaped-charge  design  calculations  on  most  computing 
systems  (large  or  small),  this  calculational  scheme  is  an  economical 
approach  for  warhead  design  and  analysis - 

V.  LIMITATIONS  OF  BASC 

The  current  version  of  BASC  requires  several  empirically  deter¬ 
mined  constants.  The  collapse  model  requires  the  determination  of  <J> 
and  k;  their  relationship  to  the  angle  of  incidence,  i;  the  constant0 
A  for  a  confined  casing;  and  finally.  Cl,  a  constant  calibrated  to  the 
105-mm,  unconfined,  shaped  charge  for  the  time  of  collapse  near  the 
apex  of  the  cone.  Penetration  requires  a  cutoff  velocity,  Umin,  and  a 
breakup  time,  Tj. 

Breakup  time,  as  used  in  the  theory  of  DiPersio,  et.  al.7,  is  a 
quantity  determined  from  a  radiograph  of  the  jet.  It  is  treated  as 
a  constant,  but  only  approximates  that  value  for  conical  liners  with 
uniform  wall  thicknesses.  Work  by  Simon21  shows  that  breakup  time  for 
a  given  geometry  remains  a  constant  over  changes  in  the  Chapman-  Jouguet 
pressure  of  186-380  kbars.  The  codes,  as  used,  requires  a  predeter¬ 
mined  value  for  T.  for  a  given  shaped-charge  geometry.  The  value  of 
T,  scales  with  cone  diameter.  All  of  these  limitations  stem  from  the 
empirical  nature  of  the  equations  of  the  code. 

Geometrical  considerations  for  conical  collapse  may  or  may  not  be 
a  limitation  to  the  BASC  code.  Research  will  have  to  be  conducted  in 
order  to  determine  this  fact.  Since  the  initial  cone  half-angle  varies 
as  a  function  of  distance  along  the  axis,  z,  and  since  each  zonal  ele¬ 
ment  has  no  interaction  with  other  elements,  all  geometrical  considera¬ 
tions  should  be  solvable. 

Current  limitations  will  be  corrected  by  future  changes  to  the 
code.  These  changes  are  explained  in  the  next  section. 


VI.  FUTURE  MODIFICATIONS 

There  are  several  tasks  intended  to  improve  BASC.  The  major  areas 
of  research  are  the  following:  (1)  modeling  the  jet  formation  for 


several  different  materials  and  geometrical  consideration,  (2)  mods.'ing 
the  viscous  effects,  as  suggested  by  Walters,22  in  a  nature  compatible 
with  BASC,  (3)  modeling  the  breakup  of  the  jet  with  reference  to  a 
maximum  strain  to  break,  (4)  modeling  the  threshold  for  a  jet-no-jet 
criteria  for  many  more  materials. 

•  As  explained  earlier,  elements  in  the  apex  region  of  a  conical 
liner  reach  the  jet  formation  region  before  these  elements  are  accel¬ 
erated  to  their  maximum  attainable  velocity.  To  account  for  this 
characteristic  of  conical,  liner  apex  flow,  which  is  so  important  in 
determining  exact  jet-tip  velocities,  transient  acceleration  is  being 
modeled  by  an  empirical  constant  determined  from  normalized,  copper 
liner  data13.  This  constant  accounts  for  the  number  of  shock  rever- 
-  berations  that  occur  in  the  copper  liner  before  it  enters  the  flow 
of  jet  formation.  We  can  improve  upon  this  constant  by  obtaining 
experimental  collapse  data  for  a  number  of  different  materials. 

Research  conducted  by  Simon  indicates  that  the  breakup  of  the 
shaped-charge  jet  may  be  related  to  a  maximum  strain,  which  is  a 
function  of  strain  rate.  These  observations  were  made  for  copper 
liners  with  many  different  explosive  fillers  ranging  from  TNT  (C-J 
pressure  186  kbars)  to  a  high  HMX  explosive  (C-J  pressure  380  kbars) , 
hut  with  only  one  charge  and  liner  geometry.  Work  is  continuing  by 
Chou,  el.  al., 23,214  to  define  the  critical  condition  for  1  ,-eakup  models 
based  on  these  results  and  will  be  added  to  the  code  where  piecewise 
strains  and  strain  rates  will  be  calculated.  We  will  continue  to 
calculate  penetration  velocity,  11,  and  will  terminate  penetration 
according  to  a  minimum  value  of  U  as  demonstrated  by  IliPersio,  cl.  al.; 
but  we  will  explore  the  use  of  Vjmi  as  the  penetration  cut-off  criteria. 
Finite-difference  codes  may  be  applied  either  directly  or  in  a  simpli¬ 
fied  form  to  generate  a  library  of  parameters  by  computational  experi¬ 
ments.  Tins  will  assist  in  the  research  of  some  of  these  critical 
parameters  utilized  in  the  BASC  code. 

As  certain  elements  of  the  shaped-charge  collapse  problem  are 
accurately  modelled,  that  section  of  the  BASC  code  will  be  modified. 

The  code  is  considered  a  "living"  code,  constantly  being  updated  but 
applied  wirMn  its  limitations  at  all  stages  of  its  development. 

27-W.  Halters,  "Influence  of  Material  Viscosity  on  the  Theory  of  ' 

Shaped-Charge  Penetration,"  ARBRL-MR- 02941,  August  '1979.  (AD  #804(4891) 

22P.C.  Chou  and  J.  Carleone ,  "Calculations  of  Shaped-Charge  Jet 

Strain,  Radius  and  Breakup  Time,"  BRLCR-246 ,  July  1976.  (AD  #B007240L) 

2l,e7.  Carleone,  P.C.  Chou,  and  R.  Ciaaarelli ,  "Shaped -Charge  Jet 
Stability  and  Penetration  Calculations,"  BRLCR- 351,  September 
1977.  (AD  VA050117) 
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VII.  SUMMARY 


BASC  is  a  simple,  well -documented,  shaped-charge  code  that  may 
be  applied  to  many  problems  to  predict  trends  and  gross  effects. 
Difficulties  in  predicting  the  jet-tip  characteristics  still  exist 
for  some  materials,  but  future  modifications  should  correct  this 
deficiency.  The  code  is  so  structured  that  it  can  *row  and  become 
more  widely  applicable  as  modeling  improvements  arc  available. 
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APPENDIX  A 

FLOWCHART  OF  COMPUTER  PROGRAM 


& 3BHBGKT3E SSi 'u E  ■.*  .H.  ^ 


CALI  PLOTS 


/  RHOJ 

RC  (I 1=a/ -  *  RSQ  (I)  X  VJO 

^  U'  V  2.  X  AKAY  X  CK 


AI=  AI/57.3 
DA  =  D/SIN  (AI) 


onnnonoonnnoonnnnnoonnnnnonoonoononooonnnonnnnnoonnono 


PROGRAM  PASC (INPUT, OUTPUT. TAPES* INPUT. TAPE b«OUTPUT) 
8PL  ANALYTICAL  SHAPED  CHANGE  (HASC)  CODE 


PROGP AHMED  BY  JOHN  7.  HAPR1SON 


BRL  ANALYTICAL  shaped  CHARGE  CODE 

BRITTEN  IN  STANDARD  FORTRAN  IV 
SEMI-tMPIRICAL  COOE  BASED  ON  THF  FORMULA- 

l/PHl»l/PHl(l*K*RHOJ*EPS/E 


UNITS  FOR  THIS  CODE  ARE  CENTAMETERS, GRAMS, AND  MjCPOSECENDS, 

THE  CODE  BILL  CALCULATE  COLLAPSE  VELOCITY  OF  THf  LINER  AND 

VELOCITIES, MASSES. AND  ENERGIES  OF  BOTH  THE  JET  AND  SLUG, 

THIS  CODE  BILL  PREDICT  PENETRATION  AND  HOLE  PROFILE  AND  6IVE 
PENETRaTION-STANDOFF  CURVES  AND  TABLES, 


LIST  OF  VARIABLES 


alrao 

CON 
DZ 
OT7 
A  I 

Z  (I) 

T7  ( 1) 
GAMA ( I ) 
BETA ( I ) 
V<1> 
TAU(I) 
PHI (I ) 
PHIP 
RPH1P 
HPHI 

cm 

VJlIl 

VN(I) 

CO 

R  V  ( I ) 

OF  J ( I ) 
OFN  ( I ) 

R  ( I) 

DA 

DHL  m 

E(I) 

DPOINT 


ALPHA  IN  RADIANS  MAIN  27 

DETONATION  PRODUCTS  CONSTANT  <K).  MAIN  2b 

DELTA  7  DISTANCE  ALONG  THf  LENGTH  OF  THF  LINER  MAIN  29 

TIME  THE  DETONATION  WAVF.  ARRIVES  AT  EACH  Z  POINT  MAIN  30 

ANGLE  Of  INCIDENCE  MAIN  31 

DISTANCE  ALONE  LINER  MAIN  32 

time  at  Each  point  on  liner  main  33 

ANGLE  BETWEEN  AXIS  AND  COLLAPSE  DIRECTION,  MAIN  3* 

COLLAPSE  ANGLE  MAIN  3b 

COLLAPSE  VELOCITY  MAIN  3r> 

TIME  FOR  EACH  ELIMENT  TO  REACH  THE  AXIS  MAIN  37 

ANGLE  OF  LINER  BENDING.  MAIN  3b 

ANGLE  OF  EXPLOSIVE  BENDING.  MAIN'  39 

1  /PHI  P  MAIN  *,0 

1 /PHI  MAIN  41 

POINT  ON  THE  AXIS  BHEPE  EACH  PARTICLE  BILL  HIT  MAIN  4? 

VELOCITY  OF  THE  JET  MAIN  43 

VELOCITY  OF  THE  SLUG  MAIN  44 

SOUND  SPEED  OF  LINER  MATERIAL (  COPPER  CCe.395  CM/MSEC)  MAIN  44 

Flow  VELOCITY  (VP),  DIRECTED  INTO  STAGNATION  POINT.  MAIN  44 

DELTA  ENERGY  IN  EACH  JET  SEGMENT  MAIN  4b 

DELTA  ENERGY  IN  EACH  SLUG  SEGMENT  MAIN  4b 

RADIUS  OF  THE  CONE  AT  EACH  DELTA  Z  DISTANCE  MAIN  47 

APARENT  DETONATION  RATE.  MAIN  *b 

DELTA  MASS  OF  THE  LINER  MAIN  49 

AMOUNT  OF  EXPLOSIVE  ABOVE  EACH  POSITION  ON  THE  LINER  MAIN  SO 

INITIALLY  IS  THE  TOTAL  HEIGHT  OF  THE  CHARGE*  THEN  BE-  MAIN  SI 

COMES  THE  INITIATION  POINT  ON  THE  Z  DIRECTION.  MAIN  52 

RADIUS  OF  CURVATURE  OF  THE  CONFINMENT.  MAIN  53 

IF  NRAO«0.-THEN  THERE  BILL  BE  A  LINEAR  THICKNESS  OF  EXPIGMAJN  54 

MAIN  55 

INPUT  PARAMETERS  MAIN  5b 


MAIN  3 
MAIN  4 
MAIN  5 
MAIN  b 
MAIN  7 
MAIN  B 
MAIN  9 
MAIN  10 
MAIN  11 
MAIN  12 
►AIN  13 
MAIN  14 
MAIN  15 
MAIN  lb 
MAIN  IT 
MAIN  lb 
MAIN  19 
MAIN  20 
MAIN  21 
MAIN  22 
MAIN  23 
MAT  N  24 
MAIN  25 
MAIN  2b 
MAIN  27 
MAIN  2b 
MAIN  29 
MAIN  30 
MAIN  31 
MAIN  32 
MAIN  33 
MAIN  34 
MAIN  3b 
MAIN  3r> 
MAIN  37 
MAIN  3E 
MAIN  39 
MAIN  40 
MAIN  4 1 
MAIN  4? 
MAIN  43 
MAIN  44 
MAIN  44A 
►  AIN  44 p 
MAIN  45 
MAIN  4b 
MAIN  47 
MAIN  *b 
MAIN  49 
MAIN  SO 
MAIN  51 
MAIN  52 
MAIN  53 
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CARU  NUMBER  1 

I JOHN  •  1.  A  PARAMETRIC  STUDY,  2,  A  NEW  CASE. 
HEAD  >  A  HEADING  CARO. 

CARD  NUMBER  2 


HA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 


57 

5b 

5* 

bO 

61 

62 

53 

54 
«S 
4b 
47 
64 


ALPHA 

EPS 

HHOJ 

RE 

H 

COF 

RHOCON 

NPLT 

NPOS 


HALF  ANGLE  OF  THE  CONE  IN  DEGREES.  MA 

THICKNESS  OF  THE  LINER, (0  IF  UNKNOWN  OR  VARIABLE)  HA 

DENSITY  OF  THE  LINER  MA 

FINAL  RADIUS  OF  THE  CONE  MA 

HEIGHT  OF  THE  CONE.  IF  H«0.»  THEN  H  W1LI  BE  COMPUTED.  MA 
CONFINEMENT  FACTOR (0  FOR  UNCONFINED,  IF  CONFINED  ■  THICKNMA 


DENSITY  OF  THE  CONFINEMENT 

NUMBER  OF  PLOTS,  O-SKIP,  1-ALL,  2-VEL.  ♦  PENETRATION 
NUMBER  OF  (R  ,2)  POSITIONS  TO  BE  READ  IN. 

IF  NPOS  IS  0  DO  tyOT  READ  IN  (R,2>  COONIDNATES. 

CARD  NUMBER  3 


RHOC 

SO 


DENSITY  OF  THE  TARGET 
STAND  OFF  DISTANCE 

IF  SOB  0.  PENATRATION  STANDOFF  CURVES  WILL  BE  PLOTTED. 
CONSTANT  FOR  DETERMINING  HOLE  VOLUME 
IF  (CK  bO.)  CK  WILL  BE  CALCULATED  AND  OMIN  WILL  THEN 
BECOME  THE  8HN  TO  BE  PF AD  IN. 

VELOCITY  MIN.  USED  IN  THE  PENETRATION  THEORY 
IF (CK  aO.)  UHIN  WILL  BE  THE  BHN, 

BREAKUP  TIME  OF  THE  JET 
IF (  T 1  ■  0.)  T1  WILL  BE  CALCULATED. 

INITALLY  IS  THE  TOTAL  HEIGHT  OF  THE  CHARGE 
RDPT  a  RADIUS  ABOVE  AXIS  WHERE  EXPLOSIVE  IS  INITIATED, 

NEXPL  ■  EQUATION  OF  STATE  NUMBER  FOP  THE  EXPLOSIVE. 


CK 


UMIN 
T 1 

OPOINT 


MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 

MA 


NtXPL 
1 
2 
3 


EXPLOSIVE 
COMP  B 
OCTOL 
LX 


DENSITY 
1.72 
1.32 
l.*4 


DETONATION  RATE  MA 


•  R 

.65 

.677* 


MA 

MA 

MA 


a  THE  NUMBER  OF  ZONES  Z-AXIS  IS  TO  SUBDIVIDED.  IF  NbO  THE  MA 


OEFAULT  VALUE  IS  Nb6<). 


DIMENSION  E  XPLO (6 ) ,  RHOHE (6) ,  DV(6> 

DIMENSION  R 1  (100) , IMAT (100) 

DIMENSION  EPSKIOO),  OAI(lt)O) 

DIMENSION  HP (300) *PL0T (3) 

DIMENSION  PSD (50) ,  PPT(50) 

DIMENSION  CONF (4) 

DIMENSION  UJK  (?C0)  ,  DMJM200) 

DIMENSION  HEAD (6) 

DIMENSION  DMLUOO),  RV(100),  DEL(IOO) 

DIMENSION  2(100),  TZ (100) »  TAU(IOO),  E(100>*  MHH100),  DPMK100I,  MA 

1  BETA  ( 100)  ,  GAMA(IOO).  R(IOO),  C(IOO).  VJ(100>,  VM100),  OMJdUtOMA 

2  .  DMN(IOO),  DlJ(lOO),  DEN ( 1 00 ) ,  RPHl(lOO),  V  ( 100)  MA 


6V 

70 

T1 

72 

73 

74 

75 

76 

77 
76 
7V 
60 
ei 
62 

63 

64 

65 


N  66 
N  67 
N 


66 

69 

«0 


N  VI 
N  92 


93 
V4 

95 

96 

97 
96 
99 

N 1  00 
N101 
N1 02 
N103 
N104 
NI05 
M  06 
MOT 
M  OH 
MOV 
M10 
Mil 
N1 12 
M  1 3 
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COMMON 

COMMON  ALRAD*  EPS  *  RHOU*  RHOC*  RF .  RPHlO.  DTZ*  SO.  CK.  OA,  h. 
1  PT,  DEGTOR 

COMMON  AMU(IOO)*  THETA(IOO).  FUOOJj  DFdOO),  TtlOC!*  DTI100) 

1  G(IOO).  P(100>.  AUCCI*  DEL  A  (10  0)  »  DVJtlOO).  OZOOT(IOO). 

2  RSQOOO).  RC(100>.  T1 .  UMJN 

COMMON  /LPtOT /  Z.  T2,  TAU.  E.  PHI*  DPMI,  BETA*  GAMA*  R.  C*  Vj 
1  OMJ  *  OMht  DEO*  OEM*  RPH 1 •  V 

CATA  MEAD(S)  /1H>/ 

DATA  CONE ( 1 ) *  CON'F<2!  / 1 OHUNCONF I  NED.  1H>/ 

data  confq).  confui  /iohconeined  *  ihv 

DATA  PlOTUJ.PDOT  <2).PL0T<3)/4HSKlP.3HALL**HVtSL  ♦  A  EN/ 

data  explo(I)  /ghcomp  «/  *  exploi2>  /ghoctol/ 

DATA  EXPLOI3)  /AH  TNT/.EXPLO ( 4 ) /fHPfcX-DA 04/ . EXPLO <5 > /4HLX 14/ 
DATA  EXPLO(f) /4H0260/ 

OATA  PHOHi(l).  RHOHE ( 2)  /1.72,  1.B2/ 

DATA  RHOHE  (  3)  /l.t>3/*  PHOHE  ( A  )  /  i  ,  64/ *  RHONE  ( 5S /I  .8  1  /  ,RHOmE  <<>  1  / 
DATA  OV(l),  DV  (?)  *DV  (3)  *DV(A)/.7RB*  ,HAfl,.fr<30*.fleO/ 

DATA  DV (5) ,DV (M/. 8749, ,*477/ 

*»•*>•*»•••#  ••*******#*•****••*#****••*#  A*  *#**#»*••*#■(**•  **■»* 

READ  in  HEADING  ONE  TIME  FOR  EACH  PARAMETRIC  STUDY 


1  CONTINUE 

READ  (5*37)  IJOHN.HEAD 
IF (EOF (6) )  9R7R5 .2 

2  APITE  (6*35)  IJOHN.HEAD 

BEGIN  LOOP  FOP  A  PARAMETRIC  STUDY 


READ  IN  INPUT  (  2  CARDS.  BOTH  WITH  7E10.3.215  FORMAT.) 


READ  IN  LINER  AND  CONFINEMENT  INPUTS 

READ  (5.36)  ALPHA. EPS. RHOJ.RF.h.COF.RMOCON.NPLT.NKS 

READ  IN  TARGET  AND  EXPLOSIVE  INPUTS 

READ  (5*36)  RHOC *  SO ,CK , UM IN . T 1 . DPOINT ,RDPT , NE XPL . N 

SH«H 

CO*. 395 

TMjj.Ci, 

JOHN  I  *  I 
TFCUT1*0, 

TECLT?*0. 

TECUT3-0. 

TFCOT  A*  0. 

TF.  COTS«0* 

TMjCUTsO. 

TmcL  )e0. 

T  Nrl- i'2*0  • 

TmcuT  3»  0 . 

TMCLTA«C. 


MA  INI  1 A 
MAIN] 15 
MAIN]  1(5 
0  *  MAIN}  IT 
MAIM  15 
.  MAIN)  1® 

.  MAIN  1 20 
■MA1N12) 

*  VN  *MA INI  22 
MA IN] 23 
MAI N1 2a 
MAIM  250 
: MAIN12RD 
MAIN127D 
MAIN129D 
MAIN12BD 
MAIM  300 
MAIM30E 
MAIM  310 
l.?<*/MAJM320 
MAIM33D 
MAlNlSfcD 
•#***MA]  N'l  35 
MAIM  3t> 
MAIN  1 37 
MA  I M  38 
MAIN139R 
MAIN139A 
MA  INI AON 
MAIMhI 
MAIN] A2 
MAIM  A3 
MAIN145 
M  A  1 N1  *»h 
MAIM  A? 
MA  IM  46 

maim-r 

MAlhlSO 

MA) N1 SIR 

MA  1M  52 

MAINJ53 

MAIMSA 

MAIM55R 

MA1N15R 

MA  IN] 56A 

MAIN 156H 

MAIMS? 

MAIM57A 

MAIM57B 

MAIN157C 

MAIN  1570 

MA1N15TE 

MAIN]  571- 

MAIMS7G 

MAIN157H 
MAIM  57 1 
M*  IN ) 57  J 
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wM*Maiaw 
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'1  1 

T**CI)T$«0. 

HAIN1STN 

•<  ' 

DVMAXOO, 

MAIN157L 

i  1 

i 

NRAD«0 

HAINISS 

j  ; 

n*pV<N£XPL> 

MAIMSO 

I 

ICOFui 

MAINlbO 

)  ( 

IF  (COF.GT.0,1  1COF-3 

MAIKI61 

CON*. 7# 

HAINJbZ 

]  . 

I  So  2 

HAIML3 

j 

IF  (W.E6.0)  N»6S 

HAINlbA 

1 

IC»R0»*0 

MAINISS 

1 

RAD*FLQAT (NfiAD) 

MAIMbfe 

5 

VJ0*0« 

MAINIbT 

1 

THLUoO, 

MAIMG7A 

) 

TNJ*C 

HAINlfcR 

( 

tej«o. 

MAlNlbS 

! 

tkf»o. 

HA1N1T0 

7MN*0. 

MAIN17I 

-i 

THL«0. 

MAIN1T2 

ii 

SUMVjoO, 

HAIN1T3 

i 

VJMAX*0, 

HAIN174 

] 

Ui)KNAX*0n 

MA1N17S 

f 

VNAS$«0.0 

HAINJ  7b 

,« 

OHASS-O. 

HAIM7T 

j 

OHJTEoO. 

HA  IN] 7R 

DO  3  J»1 t 1*00 

MAlNlTS 

2tJ>*0« 

HA  IMHO 

3  ANU(J)»0. 

HAINlbl 

C 

HA1N1G2 

C 

CONVENT  ANGLE  ALPHA  FrON 

DEGREES 

TO  RADIANS, 

MAIN] R3 

j  c 

MAINIH4 

pi*3.im5s?6S3*s7s3 

MAIN1R5 

0FGT0R*57.?«S77SSI3I 

MAIM  OR 

ALPADoALPHA/DEGTCP 

MAIM  *7 

1 

TpRFoT  AM »5*ALRA0I 

MAINIBTA 

* 

COMPUTE  THE  HEIGHT  OF  THE 

liner 

IF  NO  GIVEN  IN  INPUT. 

maimbb 

i  C 

maim*? 

IF  (H.LE.C.)  MoRF/TAN (ALRAD) 

MAJMOO 

c 

MAIMS! 

02*H/FLOAT(N) 

MAIMS? 

DT2*DZ/P 

MAIN! S3 

c 

SFT  UP  CONSTANTS  TO  BE  IN 

LATER 

eolations. 

MAINIHA 

PHOE»CON*BHOJ*EPS 

MAIN? OS 

SA«S1MALRAD) 

MAINISb 

CAoCOS(ALRaD) 

HAIN1S7 

TAoTAN(ALHAO) 

MAlNIVb 

IF  <PF,LE.0.)  PF*TA*H 

MAIMS9 

PlPS»Pl/3. 

MA1N200 

PIP  o2» 0*PI 

MAIN800A 

SlGHA»RHOJ*EPS 

MAIN201 

c 

MA1N202 

r  . 


DETONATION  POINT  IS  THE  HEIGHT  OF  THE  CHARGF  -  HEIGHT  OF  THE  CONE,MAIN?03 

MAIN?  0* 


0P0INT«DPi-  NT»H 

IF  (PPOINT.  *.o„>  DP0INT«2.«RF 
IF  (COF.GT.O.)  CONFTKuCOF 

delcof«conftk/*o. 


MAIN2C5 

MAIN20SA 

*AlN20b 

MAIN20T 

WAlNEOP 
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)".i  :«<ji 


I 
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I 

\ 

1 
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i 
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i 
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c 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 


c 

c 

c 


COF  -  CONF INEMFNT  FACTOR  (INCREASE  IN  EXPLOSIVE  TH  CKNESS  1  .  MAI N20N 

HA  1 N? 10 

COFw 1 , ♦  «GA»RH0C0N*C0NFTKw|l,9*«  l?7R/(7.pw,60  1S6*RHGU*E  PS)  *2»*0Q3*HFMA  I N21 1 
1  /(OPCINT«2,6*.S3:  MAIN21? 

UPPE  (6.2A)  ALPHA* CON. EPS. RHOJ, CONF  i  ICOK  )  . CONFTK .PHOCON. COF  .  H.RF  .MAIN?  13* 


DZ.DT2.N 

HA  1 N? 1 *W 

WRITE.  (6.2G1  FXPlOiNEXPll .D.RHOHE (NEXPLI .DPOINT ,HDPT .RHUC . SO, 

HA  I N? IS* 

P«IN*PLOTCMPlt«l) 

HAIN216* 

RPF*PF 

HAIN217 

PHI  <n«3.0/D£GTO» 

RAIN216 

INITIAL  TIME  OF  COLLAPSE  Tw*TAU, 

MAIN21P 

Tw*,l 

HA1N220 

CONTANT  FOR  TM|  ACCELEPATION  ROUTINE. 

HAIN221 

Cl»2,«CO*fPS«HHOJ/( .3Q?a.269»8.R) 

HA1N221 A 

FT£PP»,0000U1 

MAIN222 

HA  IN? 23 

SETUP  ROUTINE 

MAIN22* 
MAI N22S 

SETUP  INITIAL  POSITIONS  (R.Z) 

HA  I N??b 
MAIN227 

HI-  INSIDE  RADIUS  OF  THE  CONE.  R  -  OUTSIDE  RADIUS  CF 

THE  CONE 

.  HAIN226 

HA  I N?  ?P 

X*EPS/CA 

HA1N230 

HA1N231 

this  is  a  setup  routine  for  a  shaped  chage  liner  with 

CONSTANT 

EPSHAIN232 

HA  I N? 33 

IF  (NROS.EO.O)  NP0S-1 

MA1N23* 

NPOSwNPOS* 1 

MA1N23S 

WRITE  (6.27) 

HA1N236W 

IF  (NPOS.GE.N)  GO  TO  5 

HAi N?37 

DO  a  U»NP0S.N 

HA I  N?3*  «* 

2  ( J)wZ(J-l) *07 

HAIN23R 

THAT ( J ) *  1 

MA1N23RA 

« (0) “7 ( J) *T  A 

HAI N?*0 

R 1 td) »R ( J1 -X 

HAIN?*1 

IF  (RI(J).LE.O.)  RI  <«J )  *0  • 

HA  I Ni *2 

IF  (RI (U) .IE.O.)  IS*J*1 

MAI N2*2A 

IF(AE'S(RI  (J)-O.O)  .LE..0001  )  HI  (J)«0.0 

HAI N2*2H 

CONTINUE 

HA  I N?*3 

CONTINUE 

main?** 

WRITE  (6.31) 

MAI N2*S* 

WRITE  (6.3*)  (I. 7(1). 1*1. N) 

HA  I N2*bW 

WHITE  (6.32) 

HAIN?*7W 

WRITE  (6.3*1  (I.H(I).I«1.N) 

MAIK24BW 

WRITE  (6.33) 

HA  I N?*9K 

WPITE  (6.3*)  (I.RI (I) .1-1. N) 

HAIN260W 

HAIN2S1 

HAIN2S2 

E  0  OF  SETUP  ROUTINE 

MAIN253 

HAIN2S* 

T2(1)*DP0INT/D 

MAIN?S5 

FPSI (1 )*EPS 

MAIN2SSA 

MAIN2S6 

START  SPACE  .  TIRE  ITERATION. 

MAIN267 

HA  I N?5H 

DO  10  1*2.  IV 

HA IN259 

0(11*0. 

MAIN260 
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VN(1)»0. 

vj<n»o. 

ujk<i>«o. 

FIND  HALF  APEX  ANGLE-ALPHA 

02«7m-2(I-n 

DBuPtl j-R(j-l) 

AlPAO»ATAN2 (DRiDZ) 

SA.SIN(ALHAO) 

CAbCOS  t ALRAD) 

TAbTAN(ALRAO) 

IF ( IF  AT ( I ) ,EU  •  1 )  COb.392 
IF(lHATm.EO.l)  RHOJbb.9 
EPSb<B(I)-RI(I))*CA 

GO  TO  SUBROUTINE  CAL A I  TO  CALCULATE  INCLINATION  ANGLE 
(DA) -DETONATION  COMPONENT  TO  VELOCITY*  AND  (RPHIO)-l. 

CALL  CALAI  (I *Z  *R«0P0INT  *RO°T  * JOHN I • A I ) 

DA J (I)bDA 

COF  IS  FACTOR  FOR  CONFINEMENT. 


IF  RAO  IS  GIVEN  NO  EQ.  0*  CLACULATE  CURVED  THICKNESS  C 

IF  tPAO.GT.O.)  E ( I ) «SORT t ABS (RA0*#2- (H-Z ( 1 ) ) ••?) ) -RAD< 

CALCULATE  EXPLOSIVE  THICKNESS 
E ( I ) >COF* (RF-R (I  I )  ♦ <  COF— 1 • )*RHOCON*CONFTK/7,B 
CALCULATE  1 ./PHI 

TZ  ( I )*2 ( 1 ) / (OA*CA) 

TZ(I)bTZ(1)»TZ(1) 

THE«»Ob.S»PI-AI 
THEbALPAO-THERO 
EPSI (I)bEPS/COS«THE90I 
E<I)bE(I)/COS(THE) 

RHOEbCON*RHOU*EPSI(I) 

CONEbRHOE«C1 

RPHI (I)«RPHI0*RH0E/E(1) 

SUBROUTINE  TO  CALCULATE  the  ACCELERATION  OF  THE  LINER* 
ROUTINE  KILL  ITERATE  ON  TAU  AND  PHI, 

DO  t  Jb 1 « 30  0 

CTb£XP< (SURT<c0NE/(E(I)*TW)  )  )  ) 

RP(U)bRPHI (I)*CT 

TWbP<X)/(DA«  (SIN(ALfiAD«UO/RP(J)  )-SA)) 

IF  (J.EG.l)  60  TO  h 

IF  (ABS  (PP  (*J» -PP  (J-l ) )  .l.E.ETERM)  QO  TO  7 
t>  CONTINUE 

END  ITERATIVE  PROCESS 
7  CONTINUE 


MA] 

MAI 

MA] 

HA] 

MA] 

MA] 

MA] 

ha: 

MA] 

MA] 

ma: 

MA] 

ma: 

ma: 

(Ail,  ma: 

/phic.  ma: 

ma: 

ma: 

ma: 

ma: 

ma: 

ma: 

ma: 

of  explosive, ma: 

MA] 

♦  (H— 2 1 1 ) I  * T  A  MA] 


N2E1 
N2t>2 
N?63 
N26fe 
N26? 
N26t 
N269 
N2T0 
N271 
:N272 
,N273 
N273A 
N27* 
,N?7S 
N27G 
N277 
N276 
N279 
N2B0 
N2fal 
N2B2 
n?b  3 
N2*A 
N2BS 
N2B6 
N267 
N2BP 
N?tl9 
N290 
N291 
N292 
N293 
N294 
N29S 
N?9e 
N297 
N29t> 
N299 
N300 
N301 
N302 
N303 
N30A 
N30S 
N306 
N307 
N3  Ob 
N3C9 
N3I0 
N311 
N312 
N31S 
N3I6 
N317 
N31B 
N319 
K320 
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c 

MA1N321 

c 

The  END  OF  THE  inverse  velocity  gradient  subroutine . 

MA1N322 

c 

HAIN323 

TAU(I)»T» 

MAIN32A 

rphi  m«RP(j> 

MA1N32S 

c 

'  MAIN32N 

PHI (I)*1.0/RPHI (I) 

MAIN32T 

OPHI(I)«AeS(PHI (I-l)-PHI (I) ) 

MAIN32P 

IF  (Pi  (l>  .Lf.tl.c)  00  TO  10 

MA1N32SA 

SP«S1N(PHI  (Ill 

MAIN32V 

TPbYAN(PhI (1) ) 

MA1N330 

sap«simalrad*phi  <n> 

MAIN331 

c 

CALCULATE  TaNI&ETa-ALPhA) 

MA1N332 

c 

HA1N333 

T ALPHAS* (D7« (SAP-SA)*TPAp (1 )«DPHI t 1 ) *CA ) / EDI* ISAP-SAI -R ( 1 ) *DPHi tl )MAJN33fc 

1  *CA*TP) 

MAIN336 

c 

MAIN336 

c 

solve  fop  the  angle  beta  in  radians. 

MAIN33T 

c 

MAIN33H 

PFTA (I )*ATAN(TAlPHAB) *alrad 

MAIN33V 

IF  (I.EO.l )  SETA(l)*ALflAO*PHI(l) 

HAIK340 

c 

MAK.3A1 

c 

SOLVE  FOP  THE  POINT  OF  COLLAPSE  ON  THE  AXIS,  STAGNATION  POINT 

♦  SH.MAIN3*.? 

C  l ; ) *2 ( I ) *SP/ (CA* ISAP-5A) 1 

MA1N343 

c 

SOLVE  FOP  COLLAPSE  VELOCITY 

MA1N3LA 

V(1)«2.0*DA*S!N<PHI (S)«.S) 

MA1N34S 

GAMAU)«.5»PI-(AL«AD*.5*PHI  (D) 

MAIN346 

VJ(I)«V(I)*C0S(ALRAD*.S*PH1  (l>-.S*RETA(ln/SiN(.5*l!tTA(I)  » 

MAIN347 

VNin*VU)«SIN(ALRAD*.5*PH!  (l)».S*8ETA(n)/C0S(,S*eETA(!) ) 

MA  I  K  3**P 

r>ML(I)»PIP  *EPS»R  (I)»0Z/CA 

MAIN349 

CML(I)«DML(I)«RHOJ 

MAIN3»9A 

DLMASS«DML(1) 

MAIN3b9b 

OMN(n«.S*(UACOS(BETAim  ) 

MAIN3S0 

uMjn>«,&*(wcos(RETAU>n 

MAIN301 

AT«2.*(AlR*0a.S*PH1 (I)«.S*BETA(ln 

MA1N3B2 

0Ejm».S*(l.»C0S(AT) ) 

HAINTS3 

OEN(I)».S«a.-COS(ATn  . 

MAIN35A 

c 

MAI N3SS 

c 

COMPUTE  the  RFLATIVE  VELOCITY  OF  THE  JET  AND  SLUG, 

MAU,3Sb 

c 

THIS  IS  THE  Flo*  VELOCITY.  (VF1. 

MAIN3S7 

Rvm*.s*ivjm-vN(i>) 

MA1N3S8 

IF  (PV(I)  .GT.RVMAX)  PVMaXbRVU) 

MAIK3SBA 

c 

TML«TML*DPL<1) 

TMJ*TMO*nHJ(I)*DML(I> 

C  IF  THE  FLOW  velocity  IS  GREATER  THAN 
IF <PV (1I/CO.GT.1.23)  TMO*0. 

C  also  PEVOVE  KINETIC  ENERGY  FOR  JET. 
IF (PV (!) /C0.GT,I.23)TEJ»0. 
TMN«TMN*DMN(I)«DMLU) 

CIMJK  U)«DHJ(I)*0HU  I) 

DEL ( I > *.5*DML ( 1 ) *V ( l ) *•? 

TKE*TKE*DEL ( 1 ) 

TF.  J*TE  J+DEJ  ( I J  *D£L  1 1 1 
UJK  (I *«VJ  (II 

IF  IVJdl  .LE.VJMAXI  GO  TL  9 
VJMAXsVJ 1 1 ) 
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MAIN3S9 
MAI N3t>0 
HA1N3BL 

I »23*CO*  REMOVE  MASS  FROP  J£YMAlN3t>l  A 

MAI N3bl  B 
MA1N36IC 
MA IN3b ID 
MAIN362 
MA1N3E3 
MA  1  N3b<* 
MA  I  N3t>S 
MAINSnb 
M A I N3b7 
M  A I N  3  6  fc 
MAIN36V 


awiiuiiaiiaAtAIaiiiiSiSS  i77Hl, 


"ITT  T"*'  ’T*** 


2KT. 


1 


I 


fet 


l 


iV'V* 


1* 


ONJTf«0*JTF  ♦PMJK  i'l'i 

SU#'*  J«SUHVJ*VJ<  I)  *&**>*  <i> 

VjTt«SUHVJ/OMJTE 

V  CONTINUE 

10  eJrTorUUHE  ^  s^ct  xncpehent 


*U0N6  THE  HElfeMT  of  THE  UMfc. 


CAU  SUBHOUT XKE  TO  COMPPESS  THE  JET  »N0  ««.«»>«  HE- 
f  ii  I  VEUGB  tN»UJ**t>*'>*»ls' 

TFJCUT-0. 

DO  12  T«pim*THCUTl*DH>l*  (H 

IFIUJK(I>.6|**»>  ^r^j*tHCUT2«0HJK  (I) 

1  Ft  U JK  a  1 THCUT 3«T  HCUT  3 ♦ 0*  J*  < I > 

IF  (UJK  < ;  1 1 1  •  mcUT4*TMCUT4*DF  JK  (1  > 

IF(UJK  1  .OE..Z  Z2ruTS*T«CUT5*OHJK<n 
IFU)JK<X1  »<3E*»1 1  4«jCUT«THviCUT*DHJF.  tl > 

If  <UJK<  I  >.«..«>  Ky3,^jcUT»0EJ<:t*0ELa> 

if  wj* \\'*¥ i tout iSecun 406 j < u •«>«• <} > 

IF  |UJM1>  »6E».E>1  Itfuil«'TECUT£*PEJn)*OEWti) 

IF  lUJK  t  n  .6E  .  .* >  It  putI»TECUT 3*Df  J  u  1  •DEL  l  J  > 

IFIUJKUJ.66..3  S5S|cUT4-0EJ»l>«DEtn) 

IF  (UJKtll  *»f **?!  TFCUTSbTECUTS*DE J (1 1  *DEU  ( 1 1 

if  cujmd .SE..D  Tt:rT;!*  ™  IT 

IF  IUJKU1  .LE»UJ*jJ**l  Tr°iPUT»Ol 
tRV <11 /C0«6T#1»231  TEJCUT  0* 

irinii'/N  .«•>•»> 

,r <rv < I 1 /C0,GT ,1 *23)  TbCUTI 

5, i  ’ ;  i  .|>> 

IF  |RV  111 /C0«GT .1*231 
IF(PV<I)/CO.OT.1.23  ECUU.J 

iF<pv(n/co.6T.\.M» 

THOJT4«0. 

IHCUT3-0. 

I*CUT2«0. 

tucuti»o. 


JET 


MAIN370 
HA1N3T1 
HAIN3T2 
HA1N3T3 

HAXN3T4 

HAIN3?S> 
HAIN3TG 
HAXN3TT 
HAXNJTfe 
VEEOCITYHAIN37P 

kHM3H0 


u 


12 


13 


IF  tPV  <11 /CO .GT • 1*231 
IF  <RV til /CO. GT, 1,231 
IFIPVUI/CO.CT.  1.23> 

IF  (RVU1/C0.6T. 1,231 
IF (PV (11 /CO.OT ,1,23) 
THASS“TMASS*DHJ- i *i 
UJKRA*«UJK (II 
HN«I 

CONTINUE 

IF  <PV <11 /CO.GT.l .231 
IF  (HV  (II  /CO  «GT  »\  *231 
IF IPV (11 /C0,6T ,1 ,231 

CONTINUE 

Nl«l 

IF*^N.GE.57)  N2-5T 

CONTINUE 
#RITE  <<>*231 
-RITE  (**24) 
no  1*  1»M*N2 

phi  n  i  •*»«!<  h^e^r 
KE7Mli*BETMl!*0£6"0P 

I) -debtor 


tpAss«0PJKtx-n 

ujfnak«ujmx*ii 

HN*I*1 
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HAIN3H0 

NA1N3U1 

MAlN3»2 

HAXN3B3 

MAIN3S3A 

HAXN3b3fc 

NAlN3t>3C 

MA1N3H3U 

MAIN3B3E 

HA  IN3b3f 

HAIN3&* 

MAIN3H4A 

HAlS3b‘b 

HA1N3P4C 

HAIN3H40 

MAIN3S4E 

HA1N3»1> 

MAIN3664 
MAIN3BSR 
HAIN3febC 
MA  JN3f>S0 
HA1N3BBE 
MAIN38&F 
MAIN3PbG 

MtlN3fct>H 

MAiN3b51 

MAIN3B5J 

KAIN3B&K 

MAIN3B5L 

MAJN3B6 

PAXK38T 

MAIN386 

MAlKlSb1/ 

HAIN389* 

HAltOH^b 

MAXNSBIC 

HAIN390 

HAIN391 

HAIN392 

HAXN393 

HAIN304 

MA1N395- 

MAIN396F 

HA1N3«7 

MAXN39B 

MAXN399 

MAINaOO 
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*>’ Ui.--.it~  ..rfAuAVA: 


.? -J A<i tt: 


nPH!(I)«DPHI  ( I )  •rEGTCirf 

HA  I N4  0 1 

■1 

2(t)«Z<I)/H 

MAIN402 

I 

WRITE  (6.23)  I.2tn.TZ<I>.EPSl<n.E<I>.PHI<:>.BET*m.DPHl  (11. 

HAINA03W 

') 

1  RPH I  (I)  «V(?).k(I)  <  7  AU ( I ) <  C ( 1 ) 

MAINAOAW 

BETA  (  I )  «6ETA  ( I)  /DEC-TOR 

HAIN405 

1* 

CONTINUE 

HA  I N4  06 

i 

j 

IF  (N.LE.N2)  60  TO  IS 

HAIN407 

N1-N2.1 

HAIN40H 

1 

N2«N 

MAIN40V 

J 

GO  TO  13 

HAIN410 

1 

15 

CONTINUE 

HAINAll 

1 

tmj.o.o 

HAIN4  ]  1 A 

j 

M«1 

HAIN4  12 

1 

N?*N 

HAIN413 

J 

IF  (N.6E.57)  N2»S7 

HAIN41A 

i 

16 

CONTINUE 

HA1N-1S 

] 

WRITE  (6.23) 

HAINAlbW 

1 

WRITE  (6.25) 

HAIN4  1?W 

j 

CO  17  1»N1.N2 

HA  I N4 1 6 

1 

TKJ«TMJ.OMJK ( 1 ) 

HAIN416A 

1 

IFIfiV  (I)/CO.CT.1.23)TMJ«0. 

MAIN* 1«B 

WRITE  (6.23)  1.2(1) .Vd(I) .VN(1I .RV(l) »UKJ(I ) . DMN ( 1 ) .DEL ( 1 ) . Dhl.  ( I ) 

.MAINaIV* 

• 

IDMJK(l)  tTMj 

HA IN420W 

HA1N421 

'1 

Z(I)»Z(1)*H 

HAINA21 A 

17 

CONTINUE 

HAIN422 

j 

IF  (N.LE.N2)  60  TO  le 

HA1N423 

j 

N 1 *N2  *  1 

HAIN424 

j 

N?«N 

MAIN425 

j 

GO  TO  16 

MA1N426 

16 

continue 

MA1N427 

\ 

IF  ( ! ERROR .EG , 1 )  GO  TO  2C 

HAINA2TA 

IP (RHOC.EO.O.O)  GO  TO  20 

MAIN426 

I 

CALL  PENETRATION  SUBROUTINE  CALCULATES  JET  RADIUS. DEPTH  OF 

MA1N42V 

1 

rfnetration.  hole  radius. .and  also  penetration  STANDOFF  CURVES 

•  HA ]Na30 

MAIN431 

i 

CALL  PENTRAT  (N. 2 *R . T2 .DMJ.Bt TA • T AU. C .UJK .R V. V JO. GRASS. MN.MEAD . RSOMA I N. 32 

i 

1.PPT.NCD.ERSI.DA1) 

MA 1 na  33 

HAINA34 

i 

t 

HA  1 Na  35 

i 

DO  1R  I«1.N 

H A I NA  36 

CONVERT  VELOCITIES  TO  (MM/MICOOSEC )  FOR  PLOTTING  AND  REPORTS. 

MAIN437 

3 

ALSO  CONVERT  DISTANCES  TO  MM, 

MAIN43B 

i 

MAIN43V 

j 

BETA(I)»BETA(I)*DE6T0H 

HAINA40 

V J  ( I »  *VJ  < 1 >  *10 . 

HAINA41 

■J 

V(I)»V(I)*10. 

HAINA42 

\ 

C(D*c(i)«io. 

HA1N443 

U JK ( 1  ) « 10 ,*UJK ( I ) 

HAIN44A 

i 

2 (I)«Z(I)/H 

HAINaaAA 

IV 

CONTINUE 

MA I N 446 

i 

IF  (NPLT.EO.O)  60  TO  20 

MAIN446 

CALL  PLOTS  iN.UJK.DPOINT.RDPT.HFAD.PSO.PPT.NCD.NPLT .IS) 

MAIN447 

3 

j 

20 

CONTINUE 

MAINA4B 

L  JKMAXslIJKMAX*  1  0  • 

HAINAAV 

i 

PFIKT  SUMMARY  OF  RESULTS 

MAIN450 

MAINA61 
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i  ■ 
i 
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u 


4  \ 


WRITE  (6.39)  HEAP  „  _ 

WRJTI  (6.30)  TML.TMJ, TMN.TKE.TEu.TEJCUT.TMJCuT 
tf*MACH«RVM  AX/CO 

WRITE  (6*41)  TECUT1.TMCUU.TECUT2.TMCUT2.TECUT3.TMCI.T3* 
1TECUT4.TMCUT4.TECUT5.TMCUT5.AMMACH 

WRITE  (6* 35)  MN.2 (HN) .UJKMAX.TMASS 
NPOS-O 

IF  ( 1 JOHN.GT *  1 )  GO  TO  1 
GO  TO  2 

9H7&S  WRITE  (6,40) 

WRITE  (6*23)  l.TW.RP (Jt .RPHI (I) .CT.BETAll) 

IERPGR-1 

STOP 


C 

C 

C 


FORMATS 


MAlN«52t 

MA1N453M 

MA  IN*53* 

MA1N454W 

MA  1N454 A 

HA1N456W 

MAlMSS  A 

MA1N4S6 

MAIN457 

MAINaSPw 

MAlNxbO* 

MA1N461 

MAIM*62 

HAIN*t>3 

MAIN464 

MAIN.65 

MA/W466 


I*  FORMAT  (2Hll!s^2HEz!8xi2HTZ.ex.3HEPS*7X.2M  E.bX.SHPMl.TX.AHbETA.fcMMMfeT 
lX.4hDPHI.6X.4MRRHl.6X.AH  V.bX.AH  *  »4*  -J 3^9.7X^1*469 

25  FORMAT  (2H1I.6X.1HZ.8X.3H  VJ.bX.2HVN.7X.3H  RV.TX.3HDMJ.7X.3HDMN.7  XMAIN46V 

1 , 3HDEL.TX .  3HDML .7X.4hDMJK.9X. 3HTMJ/) 

26  FORMAT  (8E10.4)  7* 

28  FORMAT(:kH////60x!^lNPUT  PARAMETERS*//. 5X.J  ALPHA  .  •  .E6.3, 7X ,  •» K  C0MAIN67  3 

BESW.:iRa:”«iasaR  Kri,.sK-:.*eas???,:'.wsis. . «  s: 

itop  «  ,.F6  4//5X. 'LINER  HEIGHT  »  * . FI 0 .4.7X . * L 1NEP  RADIUS  •  1  .F10.4MAIN476 
W/5X.1  IHCREmIkt  OF  Z  <D2)  -  '.FlO.i.TX,'  INCREMENT  OF  TIME  (DTZ1MAIN47T 
»  it  r i n  4  TX..NUM6EP  OF  ELEMENTS  (N)  ■  • . I4.//50X. ' DETONATK N* //)  MAIN4TB 
29%0mmWwi0*m  .  '.SA.AIOOX.'DETONATION  VELOCITY(D)  .  HA1N479 

TfTn.AX.  EXPLOSIVE  DENSITY  -  .F6.3//8X.  'POINT  OF  DETONATION  FRO-  MAIN480 
2CoSi IpeltOPOXNTl .  '.FlO.I.TX.'DET.  RADIUS  (RDPT)  •  '.F1C.4//50A.MA1N461 
3 • TARGET  I N^UTS • //5X •• TARGET  DENSITY  *  •  7.3*7Xf*STAND«*0pF  DISTANCMA1N4&2 
COhItStIWI.)  -  'tF10.6//50X'MISC.  INPUTS.//SX.  MA1N463 

30*rwiIf(/J5///////.*OX.*  SUMMAPY  OF  «SULTS' /////'  LINER  MASS 

%.r\o  ftM  tthXi*JET  ►'ASS  •  ••F10#At*  G*  l«5X*|5LUG  MASS  •  *  tMAINAhfe 

?riO*4««  GM  ♦  //•  TOTAL  KINETIC  ENERGY  «  '»Flv*t»«9  EM&S#I  tOEi**/ / 9HMN4«7 

,  TfiTXi  JtT  KINETIC  ENERGY  »  '.F10.A,'  ERGS* 1 • Ot IP ' // '  TOTAL  JET  K1MA1N4G8 
4NFTIC  ENERGY  ABOVE  JET  VELOCITY  .26  CM/MSkC  ■  '.F10.4.'  ERGS"1 .OF  IMA IN4BBA 
{!"}?  tSu^jJ?  MASS  ABOVE  JET  VELOCITY  .25  CH/MSFC  -  '.F10.4,'  OMM*  N4«bH 


6  •  //) 

31  FOPMAT  <//7(3M  I.6X.2H  Z.7XI) 

32  FORMAT  <//7<3H  1.6X.2H  R.TX)) 

33  FORMAT  <//7<3H  I.6X.2HPI .7X1) 

34  FORMAT  (7(14, 4X. FT. 4.3X1) 

35  FOPMAT ('  JET  TIP  AT  I  ■  '.15,' 

1'  MEASURED  JET  TIP  VELOCITY  • 

2ASS  OF  JET  TIP  ■  ',F10.4«*  GM 

36  FOPMAT  (7E10. 3.215) 

37  FOPMAT  (1I.6A10) 

3H  FORMAT  (11.30X.6A10) 

39  FORMAT  ( 1H1 , 30X .bAl 0 ) 

40  FORMAT  ( 'ERROR'  ) 

4 1  FORMAT ( / '  KINETIC  ENERGY  ABOVE  ,6  CM/MSEC 


MAIN4BBC 
MAINAflV 
MA1N490 
MAIN4Q1 
MA.1N492 

RELATIVE  POSITION  2/H  ■  • .F10.8//MAIN493 
*.F10.B»'  MM/MICNOSEC'//'  MEASURED  MMAINA94 
*/lHl) 


1* 

2 

3' 

4 

5' 

6 

7 ' 

e 

9* 

A 

■  END 


AND  JET  MASS  ■  '.F10.4/ 

KINFTIC  ENERGY  ABOVE  •.*  CM/MSEC 


•  .FIO. 4. 

•  .F  1 0  •*  * 


AND  JET  MASS  » 
KINETIC  ENERGY 
AND  JET  MASS 


'  .FI 0.4/ 

ABOVE  CM/MSEC 

'.no.*, 


»  '  *F  1 0  *4  , 

'  KINETIC  ENERGY  ABOVE  ,t  CM/MSEC  «  '.F10.4. 
AND  JET  MASS  *  ' .FI  0.4/ 

*  KINETIC  ENERGY  ABOVE  .1  CM/MSEC  »  '.F10.4. 
AND  JET  MASS  »  '.F10.4/K 

•  MAXIMUM  RELATIVE  MACH  NUMBER  ■  '.FIO. 4/) 
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MAIN494A 

MA1N495 

MAINARb 

MAIN497 

MA1N49B 

MAIN499 

MAIN500 

MAIN501 

MA1N502 

MAIN603 
M  A  1  N  5  0  4 
MA1N505 
MA I Np  06 
MAIN607 
MAIN508 
MAIN509 
MAIN510 
MAIN611- 


1 


"***'  ^aacmgj:: :  •  TT7*.-. 


.  •  — V —  if 


*;-TiL-->y:ck  i  LMu.,  '-‘Jt-ii . 


I 


c 

c 

c 


c 

c 

c 


c 
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SUBROUTINE  PENTRAT  ( N*Z *R .TZ.DMJ* BET A *T AU *C» Vd.RV * V JO  * DMASS.Mn* 

1  HEAO.PSO.PPT.NCD.EPSI.DAI) 

COMPUTE  PENETRATION* RAO I US  OF  JET. RADIUS  OF  THE  HOLE. 

DIMENSION  HEAD (8 ) 

DIMENSION  F.PSH100).  DAI  ( 1 00 ) 

DIMENSION  PPT (60) *  PSOISO).  PVJTI50) 

DIMENSION  7000)*  TZ  O  00)  «  DP  J  (10  0)  •  BETA  (100)*  TAL(IOO).  CtlOO). 
1  VJUOO).  R  ( 1  00  >  •  HV(IOO) 


COMMON 


EPS*  PNQJi  RHOC*  RF*  PPH10,  UTZ*  SO*  CK*  DA.  H «  U* 


COMMON  ALRAD. 

I  PT.  DEOTOP 

COMMON  AMUUOO)*  THETA  (100)  «  F  ( 1 0  0 1  « 

1  G  (100) *  P(100)«  A(lO0).  DEL  A (100)  • 

Z  PSUI.00I*  RC(IOO)*  Tl,  UMIN 
PENETRATION  PHASE 
VJO»0.0 
DC  11  I*MN.N 
EPS-fPSl (I) 

0 A>DA 1(1) 

IF  II.Lt.MN)  VJ  ( 1 1  «VJ  (MM 
IF  (I.EO.l)  RETURN 
IF  (VJO.GT.O.O)  00  TO  a 
V  JO » V  J  I MN ) 

JC1«1 

AKAV*SORT (RHOJ/RHOC) 

AK 1 *AK  A** 1 » 

AK3* (3,*AKAY) / <2.*AKAY) 

AKOKJmADAY/AKI 
A  10Kb 1 .  /AK 1 
bHN»320, 

IF  iUMlN.GT.l.O)  BHNsUMIN 

IF  (CK.EQ.O.)  CK* { 2250 .*A»20*BHN  >*l.DE-5 
IF  (UMIN.GT.1.0)UMIN«(CK«I*OEE-i3SO. 1/20400. 
T1«T1« (,7S:l-.2B)/IVJ0-.2fl| 

MbMN 
A (I )bO. 

DELAIDbO. 

OF  ( I ) *0, 

F (I)«0. 

DVJI1 >bO. 

SUMMUsO. 

K«I*1 

TKK)bTZ(K)-TZ(1) 

DHASSbDMJ'K) 

IF  (VJ(K).wF.O.(')  Vj(K)s)..OE-Oe 
AMU  (K  )  bV,  ,'J(K) 

THET A (K ) » i  '  Y  .TAU(K) 

F (K) bTHETA  < K ) /AMD (K ) - C ( K ) /V  JO 
OF  tK )=F (K) 

SUMPUbaMUIK)  ••  <-AKAY)*DF  (K)'*SUMMU 
T  (K  )  :•  AMD  I K )  •*  tl  .*AKAY)*(TO»  ( 1  *  *  AK  AY )  *SU**MU ) 
DT (K ) bT IK  >-T0 

0  IK ) *  V  JO*  f  T  (K) /AMU  IK) -F  IK) ) 


DF(IOO).  T (100).  DT(IOO) . 
DVJ(IOO)*  TOOT  1100) « 


*  511*  1 
PENTT  2 
PENTT 
PENTT 
PENTT 
PENTT 
PENTT 
PENTT 
PfcNTT 
PENTT10 
PENTT 11 
PENTT 12 
PENTT 13 
PENTT 1 A 
PENTT 15 
PENTT 1 b 
PENTT  17 
PENTT In 
PENTT ) 9 
PENTT19A 
PENTT20 
PENTT21 
PENTTHH 
PENTT  .23 
PENTT2A 

PENTT2S 

PENTT26 

PENTT27 

PENTT  2b 

PENTT2V 

PENT) 30 

PENTT31 

PENTT32 

PENTT32A 

PENTT32B 

PENTT33 

PENTT3A 

PENTT35 

PENTT3b 

PENTT37 

PENTT3B 

PENTT3Y 

PENTTAO 

PENTT  A 1 

PENTT  a2 

PENTTA2A 

PENTTaBP 

PENTTA2C 

PENTTA2U 

PENTTA2E 

PENTTA2F 

PENTT  A2G 

PENTTA2H 

PENTTA21 

PENTTA2J 

PENTTA2K 

PENTTA2L 


o  o  o 


1 


A(K)»Z(K>-R(K)*COT<PFTA (K) /?* ) 

DELA (K  > sA (K) 

VJ (F  >*VJ (2) 

NCD>1 

IF  (SO.EO.C.)  NCD»25 
DO  3  J«1.NCD 

IF  (NCD.GTyl)  S0»FL0AT(J)*KF*2. 

ZO»H«SO 
TO*2C/V JO 

IF  (ZO.GE.VJO*Tl)  60  TO  I 

PT«*Kl*<VJO*Tl)**AKOKI*20*»A10K-SORT  (AkAY*AK1*UHIMT1 
1  •AKOKl«2O»**lOF)-20 

V JT« ( AK 1*  <VJO*T1)*4AKOK1*ZO**A1QN-<ZO*PT>  >/(AKAY*Tl > 

GO  TO  2 

CONTINUE 

PT«AKAY*(VJO«T1-SORT(UMIN*T1*<VJO*T1*20/AKAY) > ) 
VJT«VJ0-PT/(AK4Y*TI) 

CONTINUE 

PSO(J)«F|.OAT(J) 

P»T ( J) »PT 
PPT(J)«PPT(J)*10. 

PVJT ( J) eVJT 
CONTINUE 

IF  PLOTTING  PENATRATION-  STANDOFF  CURVE  USE  2.  CD  TO 
IF  (NCD.GT.l)  NCnI-3 
IF (NCD.GT.l)  SO«FLOAT (NCD1)<*MF»2. 

IF  (NCD.GT.l)  PT«PPT (NcDl ) •• 1 

IF  (NCD.GT.l)  VJT«PVJT(NCDI)*.l 

70aS0*H 

TO»ZO/VJO 

RETURN 

CONTINUE 

TZ(1)«TZ(I)-TZ(1> 

DPASS»DMJ(I) 

IF  (VJ(I) .LE.C.O)  VJ ( I ) ■ 1 . 0E-02 
AMU(J)«VJO/VJ(I) 

IF  (AMU(I).LT.O.)  return 
THETA(I)-TZ<I)*TAU!l> 

F (I).TMETA(1)/APU(I)-C(I)/VJ0 
DF(I)«F  (l)-F(I-l) 

SUMPU-AMU ( I ) •• (-AKAY ) <*DF  <  I  > *SUMMU 
T ( I ) ■ AHU  < I ) ♦* ( 1 . ♦ AKAY ) • ( TO* ( 1 ,*AKAY >  •SUHMU) 

OT ( I ) «T ( 1 ) -TO 

G ( I ) » V JO* (T ( I ) /AMU ( I ) -F ( I ) ) 

A(I)*2(I)-P(I)*COT(BETA(Z;/2.) 

0ELA(I)»A(I)-A(I-1) 

DVJ(I»«VJ(!)-VJ(I-1) 

D70DT ( I ) • ( T ( I ) -TZ ( 1 ) ) *DV J  < I ) /DT  Z-V J ( I ) *DEL A ( 1 1 /DTZ 
D20DT(I)»VJ(I)-( (T(I)-TZ(I) )#0VJ(I)/DTZ*DELA(I)/DTZ) 
IF  (T(I).GT.Tl)  GO  TC  6 
IF  (Z0.GE.VJ0*T1>  GO  TO  5 

CALCULATE  DEPTH  OF  PENETRATION 

P(1)«Z0*< (T(I)/TO)**AKOK1-1.0) 
RSQ(I)«2.*P(1)«EPS»DA*DNASS/DZ0DT(I) 
BSQ(I)*SOPT(ABS(PSCl(I)  )  ) 


PENTT  4?P 
PENTT42N 
PENTT*20 
PENTTA3 
PENTT** 
PENTT45 
PENTT46 
PENTT*? 
PENTT48 
PENTT  49 
«(VJ0»T1)*  PENTTSO 
PENTTS1 
PENTTS2 
PENTTS3 
PENTTS* 
PENTTbb 
PENTT5b 
PENTTS7 
PENTT5B 
PENTTS9 
PENTTbO 
PENTT61 
PENTT63 
PLOT  HOLE  PH0PENTT64 
PENTT64A 
PENTTbb 
PENTTbb 
PENTT  bT 
PENTTbb 
F'ENTTbV 
PENTT70 
PENTT71 
PENTT72 
PENTT7S 
PENTT7B 
PENTT77 
PENTT7B 
PENTT79 
PENTT79A 
PENTTBO 
PENTTB1 
PENTTB2 
PENTTB3 
PENTTB4 
PENTTBB 
PENTTBb 
PENTT67 
PENTTBB 
PENTTB9 
PENTT90 
PENTT91 
PENTT92 
PENTT 93 
PENTT94 
’  .NTT9S 
'CWTT9B 
PENTT«7 
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«.  *'  --*4;  ■*>  i  &W JU 


t  >.-.u 


CALCULATE  RADIUS  OF  THE  JET 

RC<l.)»S(jPT<RHCJ/<2.»AKAY*CK)  >4SSRT<T1/T0>*VJ0*(20/<Z0*P<1)  n«*AK3 
1  *PSO<I) 

GO  TO  9 
5  CONTINUE 

P<I>«VJO*<T(I)-TO>*Tl/<Tl*T(I)/AKAY> 

VjPeVJO* (Tl*TO/ANAY)/ (T1*T (I) /AKAY) 

DVJPeV JP-VJPO 

0200T ( 1 ) ■  (T l»T ( I ) )*DVJP-VJP*OELA(l) 

070DT(I)«-DZ0PT<I> 

PSO  1 1 ) «e.*R  ( I ) •EPS*OA«DMASS/DZODT ( I ) 

PSQ(I)*SGRT (PSO(I ) ) 

HC(I)«SORT<RHOJ/<2.*AKAY*CK) > 4RSC ( 1 ) *VJO* < 1 ,-P t 1 ) / I AKAY«V JO*T 1 ) ) 
VjPOeVJP 
GO  TO  10 
b  CONTINUE 

IF  <P(I) .LT.PT)  GO  TO  S 
IF  <JCl.Nt.ll  GO  TO  7 
JC1*2 
TP*T  <  1 ) 

7  CONTINUE 

IF  (T(I).GT.TP)  RETURN 
U  CONTINUE 

DZODT  (n»<Tl-TZ(I>)*OVj<I)/OTZ-Vj<I)*OELAU)/OTZ 
P  <1 )»<AK1« (Tl/TO) ••AK0Kl*T<I) /<T  <  1  )  ♦ AKA Y«T 1 » - 1 . 0 ) *?C 
BSO  < 1 ) «2.*R  < I ) •EPS*DA*DMASS/DZODT  < I ) 

RSO<n*SOPTtABS<RSO<IU  ) 

PC ( I) *SORT(HHOJ/<2.«AKAY*CK) > • < V JO/ AK AY > • < AK1 • (20/ I  V J0«T1 ) ) ••A10K 

i  -<zo*p(in/<vjo*Tin*Rso<n 

9  CONTINUE 
VJPOeVJ(I-l) 

10  CONTINUE 

11  CONTINUE 

SPITE  <6. 31)  HEAD 
SPITE  (6.2?) 

SPITE  16.29)  T1.TO.ZO.PT.VJT.VJO, AKAY. UN1N.CK 

M«=l 

f,?«N 

IF  (N.GE.57)  1,2*55 
SRITE  (6.30) 

12  CONTINUE 
SPITE  (6.23) 

WRITE  (6.27) 


Tl.TO.ZO.PT, VJT.V JO, AKAY.UN1N.CK 


WRITE 
ll«0 
00  13 
SPITE 


00  13  1-N1.N2 

SPITE  (6.23)  I  »  APU  (I)«THETA(I),F(I),DF<1).T(1),DT(1)»G(I)»RSC(I). 
1  A(1),0ELA(1),0VJ<I).0Z0DT<I) 

13  CONTINUE 

IF  (N.LE.N2)  GO  TO  1 A 
N 1  *N? ♦  1 
N  ?  *  N 

GO  Tp  12 
1*  CONTINUE 

00  IE  I=RN,N 
IE  CONTINUE 


PENTT9B 
PENTT99 
PENT1 00 
PENT1 0) 
PENT  102 
PENT103 
PENT1 04 
PENT105 
PENT106 
PENT! 07 
PENTlOb 
PENT  109 
PENT  110 
PENT1 11 
PENT1  12 
PENT) 13 
PENTl 14 
PENTl 15 
PENT1 lb 
PENTl 17 
PENTl  1  6 
PENTl  19 
PENTl 20 
PENT  121 
PENT122 
PENT  1 23 
PENTl 24 
PENT125 
PENTl 2b 
PENTl 27 
PENT  1 26 
PENT  1 29 
PENT  130 
PENT  1 31 
PENTl 32 
PENT  1 33W 
PENT l 34R 
PENTl 35S 
PENT136 
PENT  1 37 
PENT13P 
PENTl 39s 
PENT140 
PENT  1 4 1 S 
PENT142S 
PENTl 43 
PENT  1 44 
PENT145S 
PENT  1 46S 
PENTl 47 
PENTl 46 
PENT ) 49 
PENT  150 
PENT  1 51 
PENT  1 52 
PENT  1 53 
PENT  1 54 


missKsammssaKBESlBSi  TV, :  V..C.jii'SS0fciS3*Bkaaffivi, 


non 


Mai 

PE NT 155 

K2»K 

PENT 166 

IF  (N.GE.50)  N?a«,0 

PENT157 

(SUITE  (6*31)  HEAP 

PENT) 5BW 

S0»S0/(2.*HF1 

PENT159 

WPITE  (6*26)  SO 

PENT160W 

WRITE  (6*23) 

PENT161W 

WRITE  (6.2s.) 

PENT162W 

16 

CONTINUE 

PENT163 

IF  (M.GT.25)  WRITE  (6.26) 

PENT164W 

DO  19  I-N1.N2 

PENT165 

IF  (F(I).OT.PT)  00  TO  20 

PENT) 66 

IF  (P(l).ECl.O.C)  1 1 » 1 1  *  1 

PENT167 

IF  (I.GT.25)  GO  TO  17 

PENT166 

WRITE  (6*24)  I.P(I)  *RC(I)  *PPT  (D.PSOII) 

PENT  1 69W 

00  TO  IP 

PENTITO 

17 

CONTINUE 

PENT171 

WRITE  (6*24)  I «P ( I ) ,RC ( I ) 

PENT172W 

1H 

CONTINUE 

PENT1 73 

19 

CONTINUE 

PENTIT4 

IF  (N.LE.N2)  GO  TO  20 

PENT1T5 

N1«N2«1 

PENT176 

N2«N 

PENT1T7 

CO  TO  16 

PENT17R 

20 

CONTINUE 

PENT179 

II«I1*1 

PENT160 

DO  21  I»l» I l 

PCNTI61 

PC ( I ) "PC (II) 

PENT  182 

21 

CONTINUE 

PENT  163 

RETURN 

PENT184 
PENT1 65 

penetration  formats 

PENT166 
PENT  1 67 

22 

format (40*,* initial  conditions  for  penetration*//) 

pentimp 

23 

FORMAT  (I4.2X.1P12E10.3I 

PENT169 

24 

FORMAT  (I3.7X.2(F10,5.10X*F10.S.?OX) ) 

PENT190 

2b 

FORMAT  ( 3H  I . 12X * *P (CM) » * 15X * *PC (CM ! • .23X . • PT IMM ) *  * 

15X»*S0lCD)'/)PENT191 

26 

FORMAT  <1H].3H  ltl3X.ipi.16X.iRf) 

PENT192 

27 

FORMAT  (2h  1,6X,?HMU,6X.SHTHETA.SX,2H  F.SX.3H  DF.7X. 

2H  T.BX.2HDT .BPENT1 93 

1  X.?H  G .6X. 3HP5G.BX  » 1HA.6X.4HDELA. ,X.2HDV .6X.4HDZPT/) 

PENT  1 9* 

26 

FORMAT (10X. 'HOLE  PROFILE  SO  a  *tF6.2.3h  CD.25X. 'PENTRATION  STANPENT195 

IDOFF • / ) 

PENT195A 

29 

FORMAT ( 5X . • T1  ■  * . F10 .4 .SX. *T 0  ■  ‘.F10.4.5X,*  2ft  ■  • 

»  F 1 0 . 4 , //4X .  PENT1 96 

1*  PT  a*,Fl0.4,5X.*VJT  a*,F10.4.5X.»  VJO  -'.F10.4//4X 

. *  ANAY  a  »  PENT  1 97 

2.F10.6.5X,*  UP IN  ■  *  <F1 0 .5.5X  «  *  CK  a  *,F10,5//> 

PENT196 

30 

FORMAT ( 1H1 • •  PENETRATION  PHASE*/) 

PENT199 

31 

FORMAT  ( 1H1 .30X.6A10///) 

•  PENT200 

END 

PENT201- 

74 


c 

c 

c 

c 

c 


SUBROUTINE  C  &  L  A  I  ( I .2 «R *DPO JNT .Rt)PT , JOHNI . AI ) 


SUBROUTINE  RILL  CALCULATE  INCLINATION  ANGLE  (*!>♦  DETONATION 
COMPONENT  OF  THE  COLLAPSE  VELOCITY  U'A).  A NO  WILL  CALL  CAlPHI 
SUBROUTINE  WHICH  WILL  CALCULATE  l./PHIO  (HPH10)  GIVEN  <AI). 


CIHENSION  ZCIOO).  R  1 1 00)  _  .  „„ 

COMHON  ALHAD*  EPS*  RHOU*  RHOC»  RF *  PPHlO*  DTZ*  SO*  CK* 

1  PT.  DEGTOR 
ALPHA«ALRAD*DEGTOR 

THFTA«ATANC (R(I)-RDPT)/(DPOINT*Z (I>) > 

THETA«ThETA»DEGTOR 

AI.ro. -(ALPHA-THETA) 

IF  (A1.LT.0.)  AI«180.-<ALPHA-THETA> 
call  CALPHI  (AI) 

A I  .Al /DEGTOR 
DA-0/SINCA1) 

RFTURN 

ENP 


DA* 


H*  C. 


•  701*  2 
CAL  A  1  2 
CALAI  3 
CALAl  4 
CALAI  S 
CALAI  t> 
CALAI  7 
CALAI  B 
CALAI  R 
CALAI10 
CALAI  11 
CALAI12 
CALAl 13 
CALAll* 
CALAI1S 
CALAI 1G 
CALAI 17 
CALAl )B 
CALAI1R- 
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no  nnno 


SURPOUTINE  CALPHI  < Al > 

GJVEN  AN  INCLINATION  ANGLE  (AI)  INTERPOLATE  INTO  TABLES  AND  FINu 
AN  ANGLE  PHIO.  THEN  CALCULATE  1./PH10  (RPHIO)  AND  RETURN. 

DIMENSION  TPHI (37) t  TAI (37) 

COMMON  ALRAD,  EPS*  RHOJ,  HHOCi  PF«  RPHIO.  DTZ*  SO*  CK*  DA*  H*  0* 


1  PT.  0E6T0P 

DATA  TPHI  /C*0«  .OH*  .25,  .59*  1,02* 

1  3.30*  3.70.  4 . 1 A •  A. 51*  S.08*  5.65. 

2  7. *2*  7.56.  7.73*  7.«0*  R.04,  6,17* 

3  6.06*  7.90.  7.65*  7.75*  7.56/ 

DATA  TAI  /O.C*  .02*  .05*  .15*  .23* 


SFI-.0501 
SFPHI«.20? 

X»SFI*AI 
DO  1  J»1 *37 

IF  ( T A 1  ( J)  .6T .X  )  00  TO  2 

1  CONTINUE 

2  CONTINUE 
OXIM-TAI  (J-l) 

ROX«DXl/ (TAI (J)-TAI (J-l) ) 

VuRDX* (TPHI (J) -TPHI (J-l >  >*TPHI (J-l) 

IF  (PHIO.LE..0001)  PHIO«0. 

PHIO«V/SFPHI 

RPHI 0«1 • O/PHIO 

BPHIO»PPHIO*DEGTOR 

RFTURN 

END 


•  720*  3 
CALPI  2 
CALPI  3 
CALPI  a 
CALPI  5 
CALPI  6 
CALPI  7 
CALPI  6 

1, A3.  1.76*  2.03*  2.47*  2. 66, CALPI  90 
5.96.  6.37.  6.70*  7,01.  7, 22. CALPI 10U 
6.23*  6.26.  6.30*  6.26.  6.16.CALP1110 

CALPI 120 

46,  .60*  .72*  » 84 ,  .CALP113D 
1.92*  2.03.  2.16*  2.CALPIU0 
3.57,  3.72.  3.66*  4.CALPI15D 
CALPI 160 
CALPI  17 
CALPI18 
CALPI 19 
CALP120 
CALPI21 
CALPI22 
CALP123 
CALPI2A 
CALPI25 
CALPI26 
CALPI27 
CALP126 
CALPI29 
CALP130 
CALPI31 
CALPI32 
CALPI33 
CALP13A- 


32*  *a2« 

19A,  1.06*  1 . 1 A «  1.29,  1.44.  1.54,  1.67*  1.79, 

220.  2.41*  2.61,  2.79.  3.01*  3.14,  3.25*  3.3R, 

300,  4.16,  4.29,  4.47/ 

TaI  Scale  FACTOR  100  DEO  «  5.01  COUNTS  (SFI «b .0 1 / 1 0 0. ) 

TPHI  SCALE  FACTOR  40  DEO  ■  R.06  COUNTS  (SFPH1 «R. 08/40 . ) 
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SUBROUTINE  VELGH  < IK »UJ* « Dm JK. IS ) 

DIMENSION  UJK(?00)«  OHJK (200) 

E»0  • 

I JK  *  I K“ 1 

1  KP«) 

COMPRESSION  OF  THE  JET 
DO  ?  M«IS*1JK 

IF  (UJK (Ml  .GE .UJK  i H ♦  1 )  )  CO  TO  2 
IF  (ABS  (UJK |M)«UJK<M*1) ) .LE.l ,E-S)  GO  TO  t 
TI*PMJK (H) *UJK (M) *OMjK (M*l ) *UJK  <M*1 > 
T?*EH'MJK  <M*1 1  *  (UJK  (M)-UJK  !»♦!  )  ) 

OEH»l .O/IOMJK (HO  *DHJK (H*in 
UJK (M) *  <T 1-T2 ) *OEM 
UJK (H*1)«(T1*T?)*DEM 
KP«0 

2  CONTINUE 

IF  <KP.EO.OI  GO  TO  1 

RETURN 

END 


o  75*»  <* 
WELSH  2 
WLLGH  3 
WELSH 
VELGP  E> 
WELSH  6 
WELSH  T 
VELGP  « 
VELGP  V 
VELGHIO 
VEL6PI 1 
VELGHI2 
VELGH13 
VELGP I  * 
VELGH1S 
VELGH1G 
VELGH17 

velghih 

WELSH  IV- 


SUBROUTINE  PLOTS  (N.U JK .DPOINT .RDPT .HEAD.PSO.PPT ,NCD »NPLT * J ST ART )  •  773*  5 
DIMENSION  PSO(SO).  PPT (SO) «  S(50>  PLOTS  2 

DIMENSION  B  (5000) •  HEAD  <R)  PLOTS  3 

DIMENSION  UJK (200) »  DMJK (200)  PLOTS  4 

DIMENSION  RDU),  CV(4),  *B(4).  TC(4),  PCU),  DMNJU),  VVC(4)  PLOTS  b 

DIMENSION  ZUOO).  TZ  (100 ) »  TAU(lOO),  EdOO).  PHK100).  DPMIOOO),  PLOTS  N 
1  BETA!  100)  •  GAMA  (100).  RllOO),  C(100>.  VJ(IOU).  VMIOO),  DMU  1 1  00  >  PLOTS  7 
?  «  DMN(IOO).  DEJ(I00)«  OEN(IOO),  PPHltlOO).  V ( 1 0 0 >  PLOTS  b 

DIMENSION  F  (4 ) «  6(4)  PLOTS  B 

common  alpad.  eps,  phoj.  rhoc.  rf«  pphio.  c>tz«  so*  ck.  da.  h.  d.  plots i o 

1  PT,  OEGTOR  PLOTS  1 1 

COMMON  /LPLOT/  2.  T2.  TAU.  E.  PNI .  DPMI.  BETA.  GAME .  R*  C.  VJ.  VN.PLQTS12 
1  pMj,  OMN,  DEj.  OEN.  BPhl,  V  PLOTS  1 3 

DATA  DMNJ(l),  DMNU(2)  /lOHRELATiVE  M.  4HASS>/  PLOTS14D 

DATA  RD(l).  PD ( 2) .  RD(3)  /lOHRELATiVE  D.  IOHISTANCE  EB»  lOhOM  CONFPlOTSISD 
1  AP/  PLOTS  1 oD 

DATA  RO (4)  /3HEX>/  PLOTS17D 


OaTA  CVU),  CV  (2)  .  CV13V  /10MCOLLAPSE  V.  I0HELOCITT  JN,  1  OH  C.M/*1  CPLOTS )  fcD 
IPOS/  PL0TS19D 

DATA  CV  (4 )  /3HEC»/  PLOTS200 


DATA  CV  (4 )  /3HEC>/ 

OATA  AB(l)  .  AH  (?)  »  ABO) 
OATA  TC ( 1 ) .  TC(2).  TC  13 ) 
1)>  / 

data  pcui.  pc(2).  pco) 


/I0HANGLE  BETA.  10H  (DEGREES) •  ln>/ 


PLOTS? ID 


/10HTIME  OF  CO.  10HLLAPSE 


10HCPC  StCPLOTS22D 
PLOTS23D 


Data  PCUI.  PC(2).  PCO)  /IOHPOINT  of  C.  IOHOLLaPSE  (M,  3HM)>/  PLOTS24D 

OtTA  VVC(l).  VVC (2 ) .  V VC (3 )  /10H  V.  VJ  (MM.  lOh/MlCROSEC) .  1H>/  PLOTS25D 
DATA  T JET  /4HJET >/  PL0TS2SD 

DATA  TSLUG  /SMSLUG>/  PL0TS27D 

XPAGE*1 A  PLOTS2S 

call  PLTCCB  (XPAGEtl.RU).B(5000>>  PL0TS2R 

II«N-ISTAPT  PLOTS30 

YMAX.12,  PL0TS31 

CALL  PLOT SI  (3.0.3.0,1.0,YMAX..l,.50.rtD,VVC,lf>.10.A)  PL0TS32 

CALL  PLTCCD  (l.O.Z(ISTART) .V(ISTART).II)  PL0TS33 

CALL  PLTCCD  < 1 , 0 ,2 ( 1ST ART ) . VJ ( 1ST APT ) . 1 1 )  PL0TS3* 

CALL  PLTCCO  (4,C,Z(1START) .UJK(ISTAFT) ,11)  PL0TS35 

CALL  PLCCSP  (XS.YS.UFAC)  PL0TS36 

XX*. ?S  PLOTS37 

YY»VJ(35)*.04*.S  PL0TS3S 

YY«V(35)4.044.5  PLOTS3V 

PLOTS40 

00  EXPERMINTAL  DATA  POINTS  FOR  THE  I OSMM  UNCONFINED  SHAPED  CHANGE .PLOTS4 1 

PL0TS42 

CALL  Pl.TEXP(N)  PL0TS43 

PL0TS44 

DPOINT«DPOINT*10.  PLOTSAb 

ENCODE  (21 «?,S ( 1 )  ) DPOINT  PLOTS4S* 

WDPT»ROPT*lO.  PLDTS47 

ENCODE  (17,3,5(4)  ) RDPT  PL0TS4«* 

ALPHA*  ALPAD*DF.  GTOP  PLOTS*V 

ENCODE  (20. 4.5(B)  ) ALPHA  PLOTSbO* 

RF  *PF • 1 0 .  PL0TSS1 

ENCODE  ( 1B.S.S (12)  ) HF  PL0TS52* 

EDS>sE  PS*1  0.  PL0TSS3 

ENCODE  dR.fc.Sdt-)  )  EPS  PL0TSS4* 

YY«VMAX*.2*YS  PL0TSS4 A 

XMAX«1,0  PL11TSS4B 

XX«.S*XMAX-.2*XS«7.  PL0TS54C 


CALL  PLTEXP(N) 

DP0INT»DPOINT*10. 

ENCODE  (21 «?,S ( 1 )  ) DPOINT 
HDPT»RDPT*10. 

ENCODE  (17,3,5(4)  ) RDPT 
ALPHA*  ALP  AD»DF.GTOP 
ENCODE  (20. 4, SIB)  ) ALPHA 
RF  *PF • l 0 . 

ENCODE  ( 1B.S.S  (12)  ) HF 
EPS*EPS*l(i. 

ENCODE  dR.fc.SdC-)  !EPS 

YY«VMAX*.2*YS 

XMAX* 1.0 

XX*.S*XMAX-.2*XS«7. 


sazasraggaae 


.1  -Kv.iUat 


r*  o  o  o  r»  o 


xx*o.o 

CALL  PLTCCT (•t.HEADfl) .U..1..XX.YY) 

yv*.*s*yhax 

XX«.fr 

CALL  PLTCCT  (.)  .5(11.0. .1.. XX. YY) 

YY».bO*YHAX 

CALL  PLtCCT  <.1.S(*> .0..1..XX.YY) 

XA«.5-.1»XS»1C. 

YA»-1.1*YS 

CALL  PLTCCT  J.l.S(B) .0  .  . 1 . «XA . YA ) 

XA*.5-.1»XS*Y. 

YA.-l ,3*Y5 

Call  PLTCCT  <.1,S(12> ,0.. 1..XA.YA) 

YA«-1 ,5«YS 

CALL  PLTCCT  (  .  1  .S ( 1 b ) . 0 , . 1 . . XA« YA ) 

This  NEXT  STATEMENT  BYPASSES  SOME  RLOTTINHO  ROUTINES 
IP  (NPLT.t-T.O)  00  TO  1 

call  plots i  (3..i6..i..i60...i,?o.,RD.PC.it>,n.A> 
CALL  PLTCCD  (1 .0*2 (1ST AMT  > *C ( I  START ) >11) 

call  pltccp 

call  plot s i  (3.0.3. 0. 1 .0. 1 .0. . 1 , . 1 .rd.dmnj. is.  10,0 
xxa.35 

YYeOMN (35) -.1 

CALL  PLTCCT  ( . ? .TSLUG , 0 . . 1 . .XX •  YY > 

YY«L'MJ  ( 35 )  -  .  ] 

CALL  PLTCCT  ( . 2 , T JET , 0 . » 1 . . XX, YY ) 

CALL  PLTCCD  (1«0»2(I$7ART) »DMJ (1ST AH T). II) 

CALL  PLTCCD  C 1 . 0.2 ( ISTART ), DM). ( ISTART >. II ) 

CALL  PL0TS1  (3. 0.16. .1.0.180,..). 20. »RD,Ab,16.lO.*l 
CALL  PLTCCD  (1 *0*2 (ISTART) .HETA(ISTART).II) 

1  CONTINUE 
CALL  PLTCCP 

This  NEXT  STATMENT  part  op  PLOTTING  bypass 
DO  PENETRATION  PLOTS (  HOLE  PROFILES). 

CALL  PLOTPEN  (N.S.PSH.PFT.NCO.HEAD) 

return 

C  format  <YH  (L-H)  *  . F 6.  A  « AH  MP>) 

3  FORMAT  (5h  RD  a.FB.A.AH  MH>! 
a  format  (YH  alpha  •  « Fft .2  «5H  DE6>) 

S  FORMAT  ( 6H  RF  is  <  F  b  <  A . AH  MH>) 
h  FORMAT  (7H  EPS  *  .Ffi.S.AH  PM>) 

F.ND 


PLOTSSaD 

PlOTSSaE 

PL0TS56 

PL0TS&6 

plotsst 

plotssr 

PLOTSBY 

RL0TS6CI 

PLOTSbl 

HL0TS62 

RL0TS63 

RLOTSfeA 

PLOTSbS 

MLOTSob 

!»L0TS67 

PLOTSbR 

PLOTS6H 

PLOTSTO 

rlotsti 

PLOTST2 
PLOTST3 
MlOTSTa 
PLOTSTb 
PlOTSTE 
PL0TST7 
PLOTSTB 
PL0TS7R 
PLOTbbO 
PLOTSbl 
PL0TSH2 
PLOTSKA 
PLOTSB5 
PLOTSdb 
PLOTSbT 
PL0TSH8 
PLOTSbv 
PLOTSYO 
PLOTSY1 
PLOTSY2 
PLOTS Y3 
PLOTS YA 
PLOTSVb 
PLOTSYb 
PLOTSY7 
PL0T5YR- 
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SUBROUTINE  PLOTS!  <X0AN.  Y8AS. XNA*  . VPAX.DX .DV.PT* .t-TY  .NNCX.NNCY  ♦ 
1  NCCI 

DIPENSION  PTXU).  PTYU) 


CNY»NNCY 

CNX»NNCX 

YS»Y*AX/S. 

XS«XPAX/7. 

CAUL  PLTCCS  (XRAB.Y6AR.XHIN.YPIN.XS.YS) 


MT«.l5 

YY«-.P*YS 

XX».S*XHAX-HT*CNX*XS 

CAUL  PLTCCT  (HT.PTXU) ,0..1..XX,YY) 

XX«-.B*XS 

YY«.5*YMAX-HT*CNY*YS 

call  PLTCCT  (PT»PTY(1)«1«*0»»XX.YY) 

CALL  PLTCCA  (DX.DY.XP1N.XHAX.YM1N. YHAX.4) 
CALL  LABELA  (DX.DY.XPIN.XMAX. YHlN.YHAX, 1 . 
RETURN 
END 


0, 


1.0) 


•  071*  b 

plot i  2 
PLOTl  3 

plgti  * 

PLOTl  5 
PLOTl  t 
PLOTl  7 
PLOTl  b 
PLOTl  9 
PLOTllO 
PLOTl 11 
PLOTUZ 
PLOTl 13 
PLOTll* 
PLOT  1 IB 
PLOT lib 
PLOTl 17 
PLOTllO 
PLOT  1 1P“ 
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SlJHROL'T  IKE  PLOTPEN  (N.S.PSO.PPT .NCD.HE AD ) 

DIMENSION  TSO  (Alt  TPT  (a)  .HEAD  (ft) 

DIMENSION  S(bO).  PSO(SO).  PPT (SO ) 

DIMENSION  hOL  (? )  *  HCNdOO).  PEN  (3)  «  PENT(3),  NAD  (f  ) 

COMMON  ALRAD*  EPS*  RHOJ.  HHOC*  HF .  RPHIO.  DT2 •  SO*  C* ♦  DA,  m.  D, 
1  PT,  DESTOP 

COMMON  AMU (100) ,  THETA(IOO).  F(100),  DF(IOO).  T(100).  OT  C 1 00  > • 
l  I?  (100)  *  PtlOO).  A(IOO),  DEL  A  (100),  DVJdOOi.  D7OCT1100), 

?  RSO(IOO).  RCdOO).  Tit  UMIN 
DATA  HOL(l),  HOL ( 2 )  /lOHHOLE  PROFI.  3HLE>/ 

data  pemd.  pen<2>.  penoi  /10HRADIUS-VS-.  IOHDEMh  OF  P,  HHEN, 
)M)>/ 

Data  PENT(l).  PENT(2>.  PEnT(3>  /IOHTImE-VS-DE,  1  OhPTH  OF  PEN. 
1CM)>/ 

data  radu).  pad(2)  /ichraoius  (mm,  2h>>/ 

DATA  TSOd).  TSO  (2)  /I  PHST  ANDOFF  bH(CD)>/ 

Data  TPT(l).  TPT!2;  /IOHPENETRATIOi  7HN  (MM)>/ 
pM«P  . 

DO  1  1=1. N 

IF  (Ptn.&E.PT)  SO  TO  ? 

P(I)»-P(I! 

RCN'dt—RCtI) 

IF  (RCm.6T.RM)  PM»RC  ( I ) 

1  CONTINUE 

P (N) «-P (N) 

2  CONTINUE 
N=  I 

PMb-P (N) 

XMAXslO. 

XMIN»“XMAX 

YMIN«-PM 

IF  (PM.LT.100.)  TMIN=-100. 

IF  (PM.LT.60.)  YM1N.-60. 

YMAXxO . 

YS= (YMAX-YMIN)/P. 

XSsYS 
XRAR*3.0 
YBAR«1 5 . 0 

CALL  PLTCCS  (XBAP.TBAR.XMIN.YMIN.xb.YS) 

XAb-.?*XS*.5*12. 

Y  A* YMAX  » . 3* YS 

CALL  PLTCCT  (.20,HOL(1) .0..1..XA.YA) 

XAb-.1*XS*5. 

YArYMIN-.S*YS 

CALL  PLTCCT  (.l.RAD(l) . 0. « 1 . «XA. Y A) 

XA*-,1*XS*9. 

YA*YMIN-.75*YS 

CALL  PLTCCT  (  .  1 »S d ) . 0 . . 1 . . XA . YA) 

XAb-, 1«XS*7, 

ya.ymin-ys 

CALL  PLTCCT  ( .  1  ,S (A) .0., 1 ..XA.YA ) 

ENCODE  (lN,5.T(i)  >  SO 
X  A«»  .  1  *XS*t). 

YAbYS*) .2 

CALL  PLTCCT  l.l.TIll.O.il.i XA.YA) 

XAbXMIN-.E«XS 

YAb.E#YmIN-.1*YS».S‘«1P, 


*  6V0*  7 
PLOTN  ? 
PLOTN  3 
PLOTN  * 
PLOTN  b 
PLOTN  6 
PLOTN  T 
PLOTN  6 
PLOTN  N 
PLOTMOD 
(MPLOTNl ID 
PLOTN12D 
. (pLOTN 1 3D 
PLOTN1AD 
PLOTN 1SD 
PL0TM6D 
PLOTN 1 70 
PLOTMM 
plotnir 
PLOTN20 
PLOTN21 
PLOTN22 
PLOTN23 
PL0TN2** 
PL0TN2E 
PL0TN26 
PLOTN27 
PL0TN2H 
PL0TN2V 
PLOTN30 
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PLOTN33 

PL0TN3A 
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PL0TN*1 
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plot  na  3 
plotn<.a 

PLOTNaS 
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plotnat 

PLOTNAft 

“LOTNAS 
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PLOTNbi 

PL0TN5* 

PLOTNbb 
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CALL  PLTCCT  <.l  .PEN(2> ,  1 . .0, . XA . Y A) 

XA«-.2*XS*li>. 

YA*YS«1.4 

CALL  PLTCCT (.2,HEAD(1> .0 . • 1 . t XA. Y A ) 

0**1. 

DY«1. 

CALL  PLTCCA  (DX.DY.XPIN.XMAX.YMIN.YMAX.A) 

dy*2.*dy 

UX*5»  0 

XMINb-XHAX 

XMAXbO, 

CALL  LABEL*  <DX,DY.XA'IN.XMAXtYNJN.YMAA.-10.0,-i0.C) 

oy*o. 

xma*«-xvin 

XVJNsO. 

CALL  LABEL*  (DX.DY.XPIN.XMAx.YMIN.YPAX. 10.0*0.0) 
RC(M  *0. 

P(N)«PM 

CALL  PLTCCD  (l.O.RC(l)  .P  (1)  .M 
«CN(NJ*0. 

CALL  PLTCCD  n.O.«CMl)«P(l).N) 

XA«-.1«XS*10. 

YA«.f*«YS 

CALL  PLTCCT  (,1,S(B) . 0 . . 1 . . XA , YA ) 

XA*-. 1*XS*9. 

Y  A*  YS 

CALL  PLTCCT  (.1.5(12) .0.«1.»XA.YA> 

Ya*.6*YS 

CALL  PLTCCT  (.1.5(16) ♦ 0 . . 1 . .XA. YA ) 

IF  (NCD.6T.1)  GO  TO  3 
CALL  PLTCCP 
RETURN 
3  CONTINUE 

XMAXbPSOINCO) 

YNiXmO, 

DO  A  1*1 «NCD 

IF  (PPT(I).GT.YNAX)  YMAXoPPT(l) 

A  CONTINUE 

YPAX  *  •  0 1 *YMAX 

Y»AX*AINT(YMAX)*1.0 

ymax*ymax*ioo. 

DY*.?*YRAX 

CALL  PL0TS1  (3.0.3.0.XMAX.YMAX.1,.DY.TS0.TPT.7,B.2> 
CALL  PLCCSP  (XS.YS.UFAC) 

XA-.5*XMAX-.?*XS«7. 

X  A  *  0  . 

YA*YPAX*.2*YS 

CALL  PLTCCT (.2. HFAD(l) . 0. . 1 . . XA , Y A ! 

Y  A*-l . 1 #YS 
XA*.S>*XHAX-.1*XS*10. 

CALL  PLTCCT  («1*S(S) .0..1..XA.YA) 
XA*.S«XMAX-.1*XS*9. 

YA*-1 .3*YS 

CALL  PLTCCT  (.1.S(12),0.«1.*XA»YA) 

YA*-1,S*YS 

CALL  PLTCCT  (.1  .S(  16) ,0.» 1 . .XA.YA) 

CALL  PLTCCD  (l.O.PSC(l) .PPT (l).NCD) 

CALL  PLTCCP 
PETUPN 

?  FORMAT  <6*-  SO  *  .F6.2.AH  CD > ) 

END 
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SUBROUTINE  PLTE  XP  (N) 

DIMENSION  P£XX<60>»  REX  VJ (60 )  •  PEXVO<60).  P£XVU),  TEXP(3) 
1  PEXPX (4) 

DIMENSION  MSC!3) 

COMMON  ALRAD.  EPS*  RHOJ.  HhOCt  RF «  RPHIO*  DTZ*  SO.  CN «  DA, 
DATA  TEXP(l).  TEXPl?).  TEXP (3)  /10HEXPERIMENT,  1 OHAL  DATA 


£•013. 


DATA  (PEXX ( I ) « 1*1 *9 )  /A •  0 .  4.5.  5.0,  5.5*  6.0.  6.*.,  7.0.  7 
DATA  <PEXVJ(1),I*1,9)  /7.0I.  6.66.  6.35.  6.03.  5.62.  5.01. 
1.19.  2.25/ 

DATA  (PEXVO (I ) , 1*1.9)  /?, 050 .  2.069.  2.065.  2.0. 

1 .P7P.  1.729,  1.509/ 

DATA  (PFXPX (11,1*1.3)  /5.6b.  6.94.  6.21/ 

DATA  (PEXvm  .1*1.3)  /2.1,  1.76,  1.42/ 

M*10 

I  COUNT* I  COUNT* 1 
IF  ( I  COUNT . LE • 1 )  60  TO  1 
Wf AC  (5,5)  M,(PEXX(I),PEXVd(I).I«l.M) 

M*M*1 

PEXX (M) «.H75 
PEXVJ(M)*9.13 
GO  TO  4 

1  CONTINUE 
PEXX(10)«.e75 
PEXVU(10)«9.13 
PEXPX(4)*.9I0 
PEXV (4 ) *9 . 1 3 
F  AC*. 5*  .7606 
DO  2  1*1,9 
PfXX(I)*PEXX(I)*FAC 
PEXX(I)«PEXX(I)/P 

2  CONTINUE 
DO  3  1*1,3 

PEXPX ( 1 ) *PF  XPX ( I ) *FAC 
PEXPX  (I )*PEXPX  <I)/H 

3  CONTINUE 

CALL  PLTCCD  (2,5,PEXX<1),PEXV0(1),9) 

CALL  PLTCCD  (2,4,PEXPX(1).PEXV(1) .4) 

4  CONTINUE 
YY*9, 

XX*. 6 

CALL  PLTCCT  (.l,TEXP<l),0.»l.tXX,VY) 

YY«6. 

CALL  PLTCCT(.1.BASC(1),0,.1.,XX,YY) 

CALL  PLTCCD  (2.5.PEXX(1) .PEXVJ(l) ,M) 

RETURN 
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105-MM  SHAPED  CHAGF  SAMPLE  CASF 


HOLE  PPDFILt 


l 

P(CM) 

l 

0,00000 

? 

0.00000 

* 

0.00000 

u 

0.00000 

S 

0.00000 

0.00000 

7 

0.00000 

0.00000 

•3 

0.00000 

\  0 

0.00000 

'  \ 

0.00000 

j  ’) 

0.00000 

'  1 

0.00000 

!  4 

0.00000 

>  u 

0.00000 

t h 

0.00000 

t 

0.00000 

'  *» 

0.00000 

.  '3 

0.00000 

0 

0,00000 

l 

0.00000 

0.00000 

0.00000 

4 

0.00000 

0.00000 

0.00000 

7 

0.00000 

1 1 

o.ooooo 

o 

0.00000 

'  n 

0.00000 

*  \ 

0.00000 

0.00000 

'  < 

0.00000 

4 

o.ooooo 

0.00000 

•  ly 

O.OOOOO 

*  7 

0.00000 

' 

*.08966 

>0 

*,3*083 

0 

*.79126 

1 

5,26781 

(3,7*158 

6,2*376 

■  \ 

6.76963 

>y 

7.30876 

t* 

7.87*5* 

1 

8. *6*7* 

M 

9.08121 

'> 

9.72599 

n 

1  n ,  *  o  1 3 1 

SO  ■  3.00  CD 


PC (CM) 

o.noooo 
0.00000 
0.00000 
O.OOOOO 
0.00000 
o.ooooo 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
0.00000 
.  1*3*7 

2.09910 

l.*6*?n 

1  .*3*00 
1  .*0390 
1.37219 
1 .3395* 
1 .30S93 
1  .27163 
1.23656 
1 ,?00HO 
1 .10**0 
1  .12760 


PEN7PAT10N  STANDOFF 


PT  <PP) 

311.70730 

36*. 96680 

360.93877 

395,21636 

*00.8*566 

390,868*6 

303. *5965 

382.67361) 

37  0 .033*  0 

350,286*9 

3*8.00923 

330.90030 

326.20037 

315.63678 

305.28231 

296,12601 

*86,15*0* 

275.369M) 

265.73235 

256.26*29 

2*6.9*767 

237.775*6 

228.7*1  IS 

2 1 9,83869 

21  1,062*7 


SO (CO) 

1.00000 
2.0(1000 
3.00000 
*.00«UO 
5.00000 
6,00  (100 
7.00000 

H.oono  0 

9.00000 
10.00000 
11.00000 
12.00000 
13.00000 
1* .00000 
15.00000 
1*. 00000 
17.00000 
18,00000 
19.00000 
20.00000 
21.00000 
22.00000 
23.0000(1 
2*.0t)000 

25.noooo 
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'  S 

14.33040 

.4.3804 

16.25045 
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*7 

16.??3h4 

.8611? 

.M 

17.26609 

.8227? 

18.34664 

.80866 

•  .n 

19.48506 

.60245 

1  \ 

20.67288 

•79683 

? 

21.91220 

.78874 

■ 

23.20448 

.781 16 

../* 

24.66304 

.77303 

•.s 

25.05794 

.7643? 

27.42006 

.75408 

‘7 

- 

28.9429? 

.74494 

•  H 

30.60488 

.73845 

•  'i 

32.14110 

.72330 

TO 

33.84023 

.70078 

M 

36.69544 

.69633 

?? 

37.40375 

.6818? 
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105-MM  SHAPED  CHADS  SAMPLE  CASE 


SUMMARY  Of  HP SUCtS 


l  IMFH  HASS  *  386. 8A9?  GH 


JET  MASS  *  122.9672  C.M 


SLUG  MASS  «  262.8R20  GM 


TOTAL  kinetic  energy  » 
hital  oft  kinftic  energy 


6.3894  EHGS*1 .0E12 

6.2273  F.RGS*l.on? 


HITAL  JFT  KlNtTlC  ENERGY  APOVE  JET  VELOCITY  .2S  CM7MSEC  «> 
IOTaL  vlFT  HASS  APOVE  JET  VELOCITY  .  ?S  CM/MSEC  «  AN. 2734  GM 


5.6406  EBGS*1.0E12 


'  INF  T  1C  ENERGY  AHOVE  .5  CM/MSF.C  ■  3.R664 

k  INF  TIC  fMFPGY  AMOVE  .4  CM/MSEC  »  4,4964 

KINETIC  F.MFPGY  AHOVfc  .3  CM/MSEC  »  5.3368 

KINETIC  ENERGY  ABOVE  .2  CM/NSEC  ■  6.6294 

KINETIC  FNFPGY  AHOVE  .1  CM/MSEC  ■  6,1685 


AND  JET  MASS  « 
AND  JET  MASS  * 
AND  JET  MASS  • 
AND  JET  MASS  * 
AND  v'F.T  MASS  « 


19.2223 

27.5096 

38.0850 

53.7736 

84.8073 


AY  I HuM  RELATIVE  MACH  NUMHFR 


'FT  TIP  AT  I  * 


38  RELATIVE  POSITION  7/H  «  .37000000 


F  A SURE D  JET  TIP  VELOCITY  *  7.01100675  MM/MICROSEC 


■FASlIRED  MASS  HE  JET  TIP 


6.1633  GM 
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